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Abstract. We investigate the constraints on Supersymmetry (SUSY) arising from available precision mea- 
surements using a global fit approach. When interpreted within minimal supergravity (mSUGRA), the data 
provide significant constraints on the masses of supersymmetric particles (sparticles) , which are predicted 
to be light enough for an early discovery at the Large Hadron Collider (LHC). We provide predicted mass 
spectra including, for the first time, full uncertainty bands. The most stringent constraint is from the mea- 
surement of the anomalous magnetic moment of the muon. Using the results of these fits, we investigate to 
which precision mSUGRA and more general MSSM parameters can be measured by the LHC experiments 
with three different integrated luminosities for a parameter point which approximately lies in the region 
preferred by current data. The impact of the already available measurements on these precisions, when 
combined with LHC data, is also studied. We develop a method to treat ambiguities arising from different 
interpretations of the data within one model and provide a way to differentiate between values of different 
digital parameters of a model (e. g. sign(^i) within mSUGRA). Finally, we show how measurements at 
a linear collider with up to 1 TeV centre-of-mass energy will help to improve precision by an order of 
magnitude. 



PACS. 11.30.Pb Supersymmetry - 
of known particles 



12.60.Jv Supersymmetric models - 14.80.Ly Supersymmetric partners 



1 Introduction 

The Large Hadron Collider (LHC) will be the first col- 
lider to directly probe physics at the TeV energy scale, 
the Terascale. The LHC is supposed to provide first beam 
collisions in autumn 2009. Despite its tremendous success, 
the Standard Model (SM) of particle physics exhibits a 
number of shortcomings which - according to the belief of 
many - might be remedied by new physics showing up at 
the Terascale. One very popular extension of the SM is Su- 
persymmetry (SUSY) [1]. Among the virtues of SUSY are 
the elimination of the hierarchy problem, it can provide 
natural candidates to explain dark matter in the Universe 
and it allows for the unification of the gauge couplings at 
the scale of grand unification. Since no supersymmetric 
particles (sparticles) have been discovered to date, SUSY 
cannot be an exact symmetry of Nature at experimen- 
tally accessible energies. Unfortunately, the mechanism of 
SUSY breaking is unknown. This ignorance is efficiently 
parametrised in the Minimal Supersymmetric Standard 
Model (MSSM) [2,3] by the introduction of all possible 
soft SUSY-breaking terms into the Lagrangian with min- 
imal sparticle content in a phenomcnological way. While 
the most general MSSM Lagrangian introduces around 
100 new parameters, mild assumptions on the absence of 
flavour-non-diagonal and CP-violating terms (motivated 



by the absence of strong flavour-changing neutral currents 
and electric dipole moments of the electron and neutron) 
and on the (effective) universality of the first two genera- 
tions reduce the number of parameters to 18 (MSSM18). 
Still, it is a formidable experimental challenge to recon- 
struct 18 parameters simultaneously from future measure- 
ments. An alternative but less rigorous approach is to con- 
front specific theoretical models of SUSY breaking (which 
typically reduce the number of free parameters signifi- 
cantly) directly with data. Among the most prominent 
of such models are minimal Supergravity (mSUGRA) [4, 
5,6,7,8] and Gauge Mediated SUSY Breaking (GMSB) [9, 
10,11,12]. 

If new phenomena which arc compatible with SUSY 
are discovered at the LHC - which we assume in this 
work - one of the major challenges will be to find out 
the underlying model and to measure its parameters as 
precisely as possible. Several studies have already been 
performed to investigate the precision with which SUSY 
model parameters can be derived from measurements at 
the LHC and how much is gained by combining them with 
data from the International Linear Collider (LLC) (see 
e. g. [13,14,15,16]). So far these studies assume an accu- 
racy for the used obscrvablcs which will only be attainable 
with a fairly large integrated luminosity. Thus they reflect 
the situation in which we might be in several years from 
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now or - for studies including ILC measurements - even 
later. 

In this paper we test the compatibility of various SUSY 
models with presently available data and constrain the 
corresponding parameters. Subsequently a projection of 
the present situation to the LHC era and beyond is per- 
formed to obtain a possible time evolution of the precision 
on SUSY parameters for mSUGRA and MSSM18. The fits 
are performed using Fittino [17] version 1.5.0. The spar- 
ticle properties for a given set of Lagrangian parameters 
are calculated using SPhcno version 3.0beta [18] which is 
interfaced with Fittino via the SUSY Les Houches Ac- 
cord [19,20]. Previous work into this direction is found in 
[14],[21]-[45]. 

In this paper, presently measured "low energy" (LE) 
obscrvables are subjected to a global fit of the mSUGRA 
and GMSB model based on Markov Chain Monte Carlo 
techniques. To accomplish this, we take advantage of a 
recent compilation of up-to-date theoretical calculations 
of precision observablcs within the MSSM [44]. Also, for 
the first time, we combine future LHC measurements with 
LE observables to determine their impact in particular in 
the early phase of LHC data taking and within models 
with a large number of parameters such as the MSSM18. 

This paper is organised as follows: in Section 2 we 
define and discuss the present and future measurements 
which serve as input to the global fit. We also describe 
briefly the computer codes employed to obtain precise the- 
oretical predictions as a function of the SUSY parameters 
the data are confronted with. In Section 3, we outline in 
some detail the different methods used to estimate the 
SUSY parameters from a global \ 2 variable. The advan- 
tages and disadvantages of the two main methods, Markov 
Chains and Toy Fits with Simulated Annealing, are dis- 
cussed. We describe an approach to discriminate between 
different values for discrete parameters of the models and 
illustrate this approach for the parameter sign(/Lt) of the 
mSUGRA model. Also, a new method to deal with am- 
biguities arising from different interpretations of the data 
within the same model is discussed. In Section 4, the re- 
sults of the different fits are presented. In Section 4.1, the 
constraints on mSUGRA and GMSB parameters are de- 
rived from available measurements, including obscrvables 
from K- and £?-decays, the anomalous magnetic moment 
of the muon (g — 2) M , precision electro-weak data from 
colliders and the value of the relic density of cold dark 
matter of the Universe, J7cdm^ 2 - We also determine the 
most sensitive observables, (g — 2) M and ^?cdm^ 2 and show 
the effect of their exclusion from the fit. For the best fit 
point, we calculate the corresponding mass spectra of all 
sparticles. For the first time, the uncertainties on the pa- 
rameters are converted into error bands on the sparticle 
masses. In Section 4.2, the results from fits to LHC data 
with integrated luminosities of 1, 10, 300 fb _1 are dis- 
played for an mSUGRA model point (SPSla) which leads 
to a collider phenomenology similar to that of the best fit 
point. In Section 4.3, we LE data with future LHC data. 
We also show that by this, a stable fit of the MSSM18 can 
be achieved and the masses of most sparticles can be pre- 



dicted. Finally, in Section 4.4 we investigate how precision 
measurements of sparticles at a linear electron-positron 
collider like the ILC with up to 1 TeV of centre-of-mass 
energy will turn SUSY into precision physics. The paper 
ends with conclusions in Section 5. 



2 Measurements and Predictions 

In this section, we describe the present and future experi- 
mental data which we confront with the SUSY parameter 
space. We use three different sets of measurements in an 
incremental way. These three sets are 

1. "Low energy" observables: existing experimental data 
which have the potential to constrain the allowed 
SUSY parameter space; 

2. Simulated LHC measurements: expected SUSY mea- 
surements for the parameter set SPSla at the LHC 
experiments ATLAS and CMS for three different inte- 
grated luminosities; 

3. Simulated ILC measurements: expected SUSY mea- 
surements at the ILC running at \fs = 500 GeV and 
y/s = 1000 GeV. 

These measurements are briefly discussed in the fol- 
lowing sections. Finally, the codes used for the theoretical 
calculations are described in Section 2.4. 



2.1 Low Energy Observables 

While no direct evidence for SUSY particles has been 
found to date, these particles contribute to higher order 
corrections to measured physical observables in a well- 
defined and calculable way if SUSY is realised in Nature. 

The measurements which are exploited to obtain con- 
straints on the allowed SUSY parameter space, can be 
grouped in four classes: 

1. Rare decays of B- and K-mesons; 

2. The anomalous magnetic moment of the muon; 

3. Precision measurements and the Higgs boson mass 
limit from high energy colliders: LEP, SLC, and Tcva- 
tron; 

4. The relic density of cold dark matter in the Universe. 

For reasons of comparability, the same measured val- 
ues have been used for the fit as in [44] although some 
of them, e. g. the mass of the top quark have been up- 
dated meanwhile. The exploited measurements and their 
values are summarised in Table 1. In the next sections 
these measurements are briefly described and limitations 
on their interpretation in terms of SUSY are discussed. 



2.1.1 Rare Decays of B and K mesons 

A strong constraint on new physics can be derived from 
flavour physics experiments, especially at the i?-factories. 
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Table 1: Available measurements from B-factories, kaon experiments, LEP, SLC and the Tevatron, as well as the 
measurement of (<? — 2) M and the cold dark matter relic density. Correlations amongst the electro-weak precision 
observables as given in [46] are studied for the fit to the existing measurements. No effect of the correlations on the 
allowed parameter regions is found. 
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The reasons for that are two-fold: First, the flavour struc- 
ture of the SM is remarkably exactly realised in Na- 
ture [60] . The apparent absence of CP- violation or flavour 
changing neutral currents beyond the SM severely con- 
strains models of new physics with additional flavour mix- 
ing. In this paper, we only study flavour-diagonal SUSY 
models which by construction fulfil these constraints. Sec- 
ond, the exact knowledge of branching fractions of rare 
decays, which are helicity suppressed or occur only at loop 
level with heavy particles in the loop, strongly constrains 
also flavour-diagonal models of new physics. 

While the observables used here can be precisely mea- 
sured (within the statistical limitations of the experi- 
ment), their prediction in the SM or in SUSY is often 



accompanied with theoretical uncertainties. The assumed 
systematical uncertainties on the theoretical predictions 
are listed in Table 1 . They are added in quadrature to the 
experimental uncertainties. Amongst the most important 
constraints are the recent measurements of B s oscillations 
at the Tevatron, the branching fraction B(B — > tv) and 
the inclusive branching fraction of radiative penguin de- 
cays, B — > s-f of the B meson. 

2.1.2 Anomalous Magnetic Moment of the Muon 

Although the anomalous magnetic moment of the elec- 
tron (g — 2) e = 2a e is measured approximately a factor 
of 200 more precisely than the anomalous magnetic mo- 
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mcnt of the muon (g — 2) ;t = 2a M , the sensitivity to new 
physics of the anomalous magnetic moment of the muon 
is typically enhanced by a factor of (m^/m e ) 2 sa 43 000, 
and represents a much stronger constraint. While its mea- 
surement is undisputed, there is ongoing debate about the 
exact value of the SM prediction for (g — 2) M . The reason is 
the fact that the non-pcrturbative contribution from the 
hadronic vacuum polarisation has to be extracted from 
other experiments such as low-energy e + e~ scattering at 
BES [61] or from r lepton decays [62,63,64]. Due to these 
uncertainties, the fit in Section 2.1 is performed with and 
without using (g — 2)^ as an observable. 

2.1.3 Measurements from High Energy Colliders 

The measurements of the Z boson mass and width and of 
its couplings to left- and right-handed fcrmions in produc- 
tion and decay, the hadronic cross-section on the Z pole, 
and the W boson and top quark mass serve to constrain 
the properties of particles contributing at loop level. Due 
to their high precision and due the absence of any ambi- 
guity in the interpretation of the measurement (as c. g. in 
the case of the relic density of cold dark matter) these 
measurements represent an important input to the fit. 

As outlined in [46], there are correlations within the 
LEP and SLD asymmetry measurements in the heavy 
flavour sector, respectively, and within the Z pole observ- 
ables. The effect of these correlations on the SUSY fit 
results have been tested for the baseline fit to the mea- 
surements from Table 1 as outlined in Section 4.1.1. 

In addition and for completeness, we also use the mea- 
surements of the bottom and charm quark and tau lepton 
masses and the measurement of the strong coupling con- 
stant a s as input to the fit. 

2.1.4 Limit on the SM Higgs Boson Mass 

The exclusion of a Higgs boson with SM-like properties 
below nih = 114.4 GeV at 95% C.L. represents an im- 
portant constraint on SUSY since at leading order, the 
lightest CP-even Higgs boson has a mass below mz . Only 
due to radiative corrections, its mass can be raised up to 
at most approximately 135-140 GeV [65]. 

For a fit of the general MSSM to the existing data, the 
SM limit on the Higgs boson mass cannot be employed 
since for a given m,h the gauge and Yukawa couplings of 
the lightest CP-even SUSY Higgs boson can deviate sig- 
nificantly from their SM values. Furthermore, additional 
decay channels, e. g. h — > AA may occur. 

As shown in [55], the experimental limits obtained for 
specific parameter choices can be as low as > 90 GeV 
within the general CP-conserving and even lower for the 
CP-violating case. 

In mSUGRA, however, it has been shown that such 
deviations cannot be realised [66]. For GMSB, such a gen- 
eral analysis is not available, but it has been checked that 
the model points selected by the fits in Section 4.1 do al- 
ways maintain sin 2 (/3 — a) ps 1, where tan/3 is the ratio 



of the two Higgs vacuum expectation values and a the 
mixing angle of the two CP-even neutral Higgs bosons. 
This ensures a SM-like production of the lightest Higgs 
boson and the absence of additional decay modes such as 
h — > AA. Therefore, we can safely employ the SM limit 
on the Higgs boson mass in this study. 

In principle, the full statistical information, 
i. c. CL s+ b("i/i) on the compatibility of the search 
result with a SM Higgs boson of mass could be 
exploited and converted into a contribution to the x 2 
function of the global fit. However, due to the theoretical 
uncertainty of 3 GeV on the prediction of the Higgs boson 
mass within SUSY, the use of this information would not 
have a significant impact on the results. 

The recent exclusion of the SM Higgs in a small mass 
region around 160 GeV by the Tevatron experiments [67] 
is not considered since the SUSY models under study do 
not allow mh above approximately 135 GeV [65]. 



2.1.5 Cold Dark Matter Relic Density 

The results from the WMAP satellite on temperature fluc- 
tuations of the cosmic microwave background together 
with various other cosmological constraints have estab- 
lished a cosmological standard model, in which approxi- 
mately 23% of the total energy of the Universe is con- 
tained in cold dark matter (CDM). This is expressed in 
the fit in terms of J7cdm^ 2 - While the presence of dark 
matter is relatively undisputed, its nature is still unknown. 
If SUSY is i?-parity conserving and the lightest SUSY par- 
ticle (LSP) is neutral and sufficiently heavy to contribute 
to cold dark matter (or if a metastable neutral sparticle 
exists with a lifetime comparable to the lifetime of the 
Universe), it contributes to dark matter through its relic 
density, and it can make up all or part of the cosmolog- 
ically observed dark matter. Therefore, in the fits shown 
later (see Section 2.1), the observable ^cdm/i 2 is used in 
different ways or not at all as a constraint. 



2.2 LHC Observables 

As a case study, we assume that SUSY is realised with 
parameters as specified in the SPSla parameter set of [68]. 
In Section 4.1), it will be shown that this bulk region 
point leads to a collider phenomenology rather similar to 
the best fit point obtained from low energy measurements. 

For the SPSla point, SUSY particles will be copiously 
produced at the LHC and a rather rich set of independent 
observables related to the masses and branching fractions 
of SUSY particles can be reconstructed. Many detailed 
experimental studies for this point (or phcnomcnologically 
similar points) exist. 

A direct reconstruction of SUSY particle masses at 
the LHC is difficult due to the escaping LSPs. Therefore, 
where ever possible, we use observables which can be di- 
rectly measured as input to the global fit. Such observ- 
ables are the positions of kinematic edges and endpoint 
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of invariant mass spectra. Where mass peaks can be re- 
constructed, those are used as well. Also, two ratios of 
branching fractions are employed. 

Measured production rates are not considered in this 
study for two reasons: First, the prediction of rates has 
rather large theoretical uncertainties, in particular if the 
production mechanism involves the strong force. Second, 
the calculation of the theoretical prediction - if realistic 
experimental cuts are taken into account - is very time 
consuming since usually Monte Carlo techniques have to 
be used to obtain these predictions. Furthermore, the in- 
herent statistical fluctuations of Monte Carlo predictions 
easily cause oscillations during the x 2 minimisation which 
destabilise the result. 

Three different integrated luminosities are considered 
separately to define the sets of accessible observables and 
their statistical and systematic errors: 1 fkT 1 , 10 flo" 1 
and 300 fb _1 . A centre-of-mass energy of 14 TeV is as- 
sumed throughout. Most of the statistical uncertainties 
are taken from [13]. Where ever results for the specified 
integrated luminosities are not available, reasonable inter- 
polations/extrapolations are used. Dominant experimen- 
tal systematic errors are expected to arise from the uncer- 
tainty of the lepton energy scale (LES) and the jet energy 
scale (JES). The LES uncertainty is assumed to be 0.2% 
for an integrated luminosity of 1 fb~ and 0.1% for higher 
integrated luminosity. For the JES uncertainty, 5% (1%) 
arc assumed for 1 ffcT 1 (> 1 fb _1 ). We assume that the 
energy scale uncertainties directly translate into equally 
large relative uncertainties on the positions of endpoints 
in mass spectra in case of fully leptonic or fully hadronic 
final states. Following [69], half of the relative JES un- 
certainty is assumed as uncertainty on the endpoint for 
invariant mass spectra involving both leptons and jets. 
Uncertainties on the endpoints related to the JES and 
the LES are considered 100 % correlated between differ- 
ent measurements. Table 2 summarises all employed LHC 
observables together with their assumed uncertainties. 

For the SPSla point, it is possible to reconstruct suffi- 
ciently long decay chains of subsequent two-body decays, 
such that mass information can be extracted in a model 
independent way. However, it is necessary to assign the 
observed decay products to the correct SUSY particles 
from which they originate. The full combinatorics for de- 
cay chain ambiguities is not yet considered in this analy- 
sis. We assume that all decay chains are correctly identi- 
fied. As a new approach to check the possible impact of 
misidentifications, we study the impact of a wrong assign- 
ment of an endpoint to its SUSY particles on the global fit 
as a case study. A more comprehensive analysis of these 
effects remains to be done. 

Most information on SUSY particle masses within the 
SPSla point can be obtained from the decay chain 

92 - qx°2 - q^ei - qt+rx !, (1) 

where £ denotes either electrons or muons. In total there 
are five different measurable invariant mass combinations 
possible for this decay chain: m^ ax , to^ x , to<w£, w ) 
and ™]whigh) (f° r their definition see e. g. [69]). 



Similar to m^", it is also possible to measure the end- 
point m TT where the electrons/muons are replaced by tau 
leptons. The ratio of the total number of events in the m« 
and the m TT distributions (corrected for efficiency differ- 
ences) provides a measurement of 

B(x°2 - W) x Wr -> x°l) 

B(X° 2 - ur) x B(fi - X?r) ' 1 ' 

Apart from the mu endpoint there is also one addi- 
tional mu measurement included. It originates from x° 
decays instead of x*2- Decays of X3 do not provide a visi- 
ble mu endpoint since it is mostly Higgsino for the con- 
sidered benchmark point and therefore the couplings are 
too small. There are several Xa decay chains providing two 
oppositely charged leptons. Of all possibilities, the chain is 
chosen which provides the largest mu endpoint within the 
SPSla scenario. This is the case for xl -> ~> £ + Z~x\- 
The other endpoints are unlikely to be measurable due to 
the superimposed spectra from the other di-lepton decay 
channels. 

Similar to the t — > r replacement it is also possible 
to exchange light-flavoured jets q with b-flavoured jets b. 
This yields a separate mf^l° s measurement. 

Jets which carry b-flavour also play an important role 
in obtaining information about the gluino mass. In SPSla, 
the gluino decays via g — > qq where q can be any quark 
flavour. Due to combinatorial background, gluinos can be 
reconstructed best, if one focuses on 

g -> 661,2 -> bbx° 2 -> 66^^ -> bb£ + rxi (3) 

From this decay chain, the gluino mass can be recon- 
structed by calculating the invariant mass of the x^bb sys- 
tem, provided that the X2 momentum is known. Due to 
the invisible x\i the X2 momentum cannot be measured 
directly, but in the chosen scenario it can be approximated 
reasonably well by [13] 

P(X°) * (l " —) P?t- (4) 

It turns out that the gluino mass estimate from this ap- 
proach is highly correlated with the assumed Xi mass, 
such that effectively — m^o is measured. Similarly 
Trig — m~ bi and mg — mg 2 can be determined by measuring 

the difference between the invariant mass of the x\bb and 
the X2P system. The ratio of the total number of events in 
the 61 and 62 mass peaks can be used to determine 

B(g -> b 2 b) x B(b 2 -> xlb) 

B(s->&i6)x*(&i 

The "stransverse mass" mT 2 [70,71] is used to extract 
information on the qa and the £l mass. The exploited 
decay chains are qn — ► qx\ and £l — > £x\, respectively. 
£l is studied in direct electro-weak di-slepton produc- 
tion via an s-channel Z/7 exchange. It turns out that the 
endpoint of the tot2 spectrum depends on the assumed 
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Xi mass in such a way that roughly y m | fl — ^ TO |f> an< ^ 
— 2m| are measured, respectively. 

Stop and sbottom sector information is obtained by a 
measurement of the endpoint of the invariant mass spec- 
trum of the tb system from the decay chains 

g - th - tbxf (6) 
g -> fo&i -> ibxf . (7) 

The variable to^ used in our fits is a branching fraction 
weighted average mtb endpoint for decay (6) and (7) to ac- 
count for the possibility that the two endpoints might be 
too close to each other to be experimentally distinguish- 
able. 

To reconstruct charginos, the decay chain 

q L -> qxt -> 0Wx\ -> qqqXi (8) 

is exploited. The chargino mass is obtained from the in- 
variant mass of the qqxi system where the two quarks 
come from the W decay. The momentum of Xi is recon- 
structed (up to a two-fold ambiguity) using a technique 
described in detail in [72]. 

The most precise determination of the Higgs boson 
mass for the considered mass range is obtained from mea- 
surements of the invariant mass of the four-lepton system 
in the decay h -> ZZ -> £+£-£+'£-' and the di-photon 
mass of the decay h — > 77. The top mass is measured 
from a combination of several different final states and 
techniques, the most precise of which is a kinematic fit for 
the semi-leptonic final state, where one W decays hadron- 
ically and the other W leptonically. 

2.3 ILC Observables 

At a future linear electron positron collider like the ILC, 
a huge variety of precise measurements of SUSY parti- 
cle properties from their electro- weak pair-production pro- 
cesses. 

In this paper, in order to illustrate the potential of a 
linear collider, a subset of the observables used in [15] 
is used. All expected mass measurements of [13] are 
used together with the expected measurements of ab- 
solute Higgs branching fractions and a large variety of 
cross-sections times branching fraction measurements of 
all kinematically and statistically accessible SUSY final 
states. In contrast to [15], only measurements at y/s = 
500 and 1000 GcV and at polarisations (P e - , P e + ) = 
(±80%,=f60%) are used, assuming a long running time 
of the ILC with C mt = 500 fb -1 on each polarisation at 
y/s = 500 GeV and at y/s = 1 TeV, respectively. The crite- 
ria used for the selection of expected cross-sections times 
branching fraction measurements is outlined in [15]. 

2.4 Theoretical Predictions 

Different theoretical codes have been used for the predic- 
tion of the observables. The low energy observables arc 



calculated by a selection of codes combined in the so- 
called Mastercode [44] . The RGE running of the parame- 
ters of the high-scale models down to the SUSY breaking 
scale are accomplished with SoftSUSY [73] . Subsequently, 
the observables of the Higgs sector and for (g — 2) M arc 
accomplished with FeynHiggs [74,65,56]. The flavour ob- 
servables are calculated with Supcrlso and other codes 
based on [75,76]. The electro- weak precision observables 
are derived in [77,78] and the cold dark matter relic den- 
sity is calculated by Micromcgas [79]. 

The SUSY mass spectrum for the LHC measurements, 
all direct sparticle decay branching fractions and the cross- 
sections for ILC are calculated with SPheno [18]. 

The known systematic uncertainties for the presently 
available observables are included and listed in in Table 1 . 

Systematic uncertainties for the LHC predictions are 
estimated from the difference of the predictions between 
different RGE codes and from scale variations. The differ- 
ences between RGE codes like SoftSUSY and SPheno are 
generally within the statistical and systematical measure- 
ments of the LHC measurements for C mt = 1 and 10fb~ , 
effectively making the LHC fits relatively robust against 
theoretical uncertainties on the order of 10 — 20 GeV [80, 
81], which is in the order of magnitude of the jet energy 
scale uncertainties. For £ int = 300 fb -1 the experimental 
systematic uncertainties are expected to be smaller than 
the estimate of the theoretical uncertainties, especially in 
case of rrih, however improvements on the precision of the 
predictions can be expected until the LHC has acquired 
£ int = 300 fb" 1 . Theoretical uncertainties will be included 
for all luminosities into the fit at a later stage. Neverthe- 
less we cross-checked the influence of an additional 3 GcV 
uncertainty on rrih due to unknown higher-order correc- 
tions for some of our results and found that it does not 
have a significant effect on the fit results for high-scale 
models. 

A special case is the prediction of ficnwih 2 in GMSB 
models. Since the gravitino is typically a very light LSP in 
GMSB with a mass in the order of several MeV, it repre- 
sents more hot than cold dark matter. Therefore ftcnwih 2 
is not used as an observable for GMSB. 

For the fits with ILC, theoretical uncertainties could 
play a major role, because the experimental precision as- 
sumed in [15] is smaller than the current theoretical un- 
certainties even in the gaugino and squark sector, where 
uncertainties of around 1 GeV are expected [80]. How- 
ever, the possible increase in theoretical precision until 
the existence of the ILC is yet unknown, hence theoretical 
uncertainties will be introduced into the fits with ILC at 
a later time. 

Fittino and the calculator programs for the predictions 
are interfaced using the SUSY Lcs Houches Accord [19]. 



3 Parameter Estimation 

In order to asses the consistency of a theoretical prediction 
(defined by a set of parameters within a specific SUSY 
model) for a given set of measurements the following x 2 
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Table 2: LHC observables which serve as input to the fits. The shown nominal SPSla values have been calculated 
with SPhcno. Most of the statistical uncertainties are taken from [13]. Where numbers for the specified luminosities 
are not available, some interpolations/extrapolations are used. Uncertainties on the endpoints related to the jet energy 
scale (JES) and the lepton energy scale (LES) are considered 100 % correlated among different measurements. 
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is used: 

X 2 = (M - 0(P)fcovj}(M - O(P)) + limits. (9) 

Here M is a vector containing the list of measurements, 
0(P) a vector with the theoretical predictions for these 
observables for a given point in parameter space P. covm 
is the covariance matrix specifying the uncertainties and 
correlations of the measurements M. In addition to the 
actual measurements M, limits on observables can be 
specified (e. g. the limit on the SM Higgs mass in case 
of the fit of a model which ensures the presence of a SM- 
like Higgs boson). This is incorporated for m lower (LL) 
or upper (UL) limits ^l/ll ^ e following way 

m f (0<(P) - L^ fja} for % {P) > 
limits = { ( L l L ~ Oi(P)) 2 /<Tf for Oi(P) < L\ h 
i=i { for Lf L > Oi(P) > Ll L 

where <Tj specifics how steeply the limit is rising once it is 
reached. 

Being measurements M and cov^ 1 are independent of 
the theoretical model they are confronted with. Contrary 
to that 0(P) depends on the model and it even depends 
on the interpretation of the data within a given model 
due to ambiguities in the mapping of an observed final 
state to its physical origin within the model (e. g. the 



assignment of kinematic edges in LHC mass spectra to the 
decays of the respective SUSY particles). The covariance 
matrix covj\/ is the sum of the statistical, systematical 
and theoretical covariance matrices, where the former is 
diagonal for independent measurements and the latter two 
can contain off-diagonal elements describing correlations. 
Using the \ 2 expression of Equation 9, the following tasks 
can be addressed: 

— find the absolute minimum i. e. the parameter 
point of a given model which fits the data best; 

— determine the P-value of the data given a best fit pa- 
rameter point of a model; 

— find secondary minima which could be confused with 
the absolute minimum; 

— derive the probability that a secondary minimum of 
the x 2 surface of the exact observables in a given model 
turns into the absolute minimum of the experimentally 
observed x 2 surface due to statistical and systematical 
uncertainties of the experimental observables; 

— derive the parameter uncertainties and correlations 
around the absolute minimum, with and without tak- 
ing ambiguities in the interpretation of the data into 
account; 

— derive the probability that due to the statistical and 
systematical uncertainties of the experimental data the 
true model of new physics is yielding a worse V- value 
than an alternative, wrong model of new physics; 
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— derive predictions for most probable values of observ- 
ables (and their expected variations) which are not 
used in the fit. 

In the following, the statistical techniques used for these 
tasks are introduced. Their application is described and 
their advantages and shortcomings are discussed. We also 
propose an approach to estimate the uncertainties on pa- 
rameters in the presence of different ambiguous interpre- 
tations of the data within the same model. 

In this paper, Gaussian uncertainties are assumed both 
for statistical and systematic uncertainties. For system- 
atic and theoretical uncertainties, there are other possi- 
ble choices. For example, box-shaped contributions to the 
\ 2 (instead of a quadratic function) are investigated else- 
where [45]. Given the general uncertainty on systematic 
and theoretical errors, we assume here that the final re- 
sult does not depend on such subtle differences. In fact, a 
larger effect can arise from unknown correlations among 
the systematic and theoretical errors. 

The x 2 hyper-surface for all considered SUSY models 
is highly non-trivial. As already shown in [17], gradient- 
based algorithms for global minimisation like MINUIT [82] 
are insufficient for most of the problems under study. 
Rather more elaborate methods, based on Markov Chain 
Monte Carlo and Toy Fits algorithms are exploited to ef- 
ficiently scan the multi-dimensional parameter space. 



3.1 Minimisation and Scanning Techniques 

Two different parameter estimation techniques are used 
in the following, a Markov Chain Monte Carlo and Toy 
Fits. These are briefly described in the following sections. 
Strong emphasis is laid on ensuring that the global mini- 
mum is found, that the errors are accurate, that the result 
is stable against different starting values, and that the 
sampling of the parameter space is fine-grained enough. 
This means that as many A^-dimensional parameter com- 
binations (where N is the number of parameters of the 
problem) as possible are actually scanned at least within 
the range of x 2 — Xmin < 6 (approximately corresponding 
to the two-dimensional 95 % uncertainty interval around 
the best fit point). These two techniques are chosen be- 
cause they are complementary in the way the uncertainties 
arc defined and in the assumptions made for the definition 
of the uncertainties. An agreement in the uncertainties be- 
tween the two methods provides a further strong evidence 
for the validity of the result. 



3.1.1 Markov Chain Monte Carlo 

The advantage of the Markov Chain Monte Carlo method 
is that it allows to obtain an efficient scan of the x 2 surface 
around its minima. Furthermore, it can be easily arranged 
that the sampling density in parameter space directly pro- 
vides a likelihood distribution for the SUSY parameters in 
the Bayesian approach (see e. g. [39], [41]). 



A Markov chain is a sequence of points Xi (i = 1, . . . , n) 
in parameter space. Each of these points Xi has an asso- 
ciated likelihood £(xj). For our study we use 

C = exp(-^y (11) 

Using the Metropolis algorithm [83], a new point x n +i 
which is randomly chosen according to a proposal proba- 
bility density is added to the chain if C(x n +i) > C(x n ). 
Otherwise it is accepted with probability C(x n +i) / C(x n ) ■ 
If the new point x n+ \ is not accepted, the old point x n 
is added to the end of the chain again and the process 
continues. The result is - under weak assumptions - inde- 
pendent of the specific choice of the proposal probability 
density function in the limit of infinite statistics. However, 
for finite statistics (even for order of 10 million parame- 
ter points for a typical 9-paramctcr model) the efficiency 
of the sampling strongly depends on the proposal distri- 
bution. Fittino implements the choice of box-shaped or 
Gaussian proposal distributions, where the width of the 
box or Gaussian can be adapted for each parameter. For 
each model and observable set, a set of pre-runs with sev- 
eral thousand points per chain is used to adapt the width 
parameters of each parameter individually such that the 
ratio of accepted and rejected points in the chain lies be- 
tween 0.8 and 1.2, for which the best scanning efficiency is 
expected. This procedure takes the initial uncertainties on 
each parameter from the pre-run into account and is re- 
peated manually until the result converges. For the results 
presented in this paper, only Gaussian proposal distribu- 
tions are used. 

The resulting Markov Chain can be interpreted in two 
different ways. In the Bayesian interpretation, it can be 
shown [84] that, if the proposal probability density is prop- 
erly chosen, the sampling density of points x% is propor- 
tional to the likelihood distribution £, which in turn is 
proportional to the posterior probability in the case of flat 
priors (as assumed in this paper). Therefore, the best fit is 
obtained at the parameter point with the highest sampling 
density oc £ max . The error on an individual parameter 
(or a subset of D parameters) is derived by integrating 
("marginalising") the sampling density over all parame- 
ters apart from the parameter(s) under study. The result- 
ing /^-dimensional distribution can then be interpreted in 
terms of — 2 ln£ + 2 ln£ max , where £ max is the likelihood 
for the parameter point with the highest likelihood. The 
1 a uncertainty of a one-dimensional parameter distribu- 
tion is defined by the region within —2 In C+2 ln£ max = 1. 
In this interpretation, the marginalised C is the proba- 
bility distribution of the true parameter value given the 
measurement. 68 % of this distribution is contained within 
la. 

The Bayesian interpretation has to be handled with 
care for two independent reasons. First, the outcome can 
have a strong dependence on the chosen prior probability. 
This is c. g. exemplified in [85]. Second, for very complex 
parameter spaces with many parameters (typically 8 to 18 
in the case of the fits presented here) one needs to check 
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Fig. 1: tan/? sampling behaviour of Markov chain using LE and LHC measurements for an integrated luminosity of 
1 fb _1 within mSUGRA for two different starting points. 



carefully that the sampling has not only reached all rel- 
evant areas of the parameter space, but that in addition 
the sampling is completely in equilibrium, i. c. that the 
likelihood is really proportional to the sampling density. 
This problem is exemplified in Fig. 1, which shows the ini- 
tial behaviour of two Markov Chains scanning same model 
space (mSUGRA) and the same measurements but with 
different starting points in the parameter tan [3. While 
above n f=a 3000 no dependence on the starting point can 
be observed in this example, below n ss 3000 the point 
density is obviously not proportional to the likelihood. 
Therefore including this region into the calculation of the 
point density would distort the result unless the Markov 
Chain length is so large that this region has negligible im- 
pact. For all results presented in this paper which employ 
the Bayesian interpretation of the Markov Chain, two dif- 
ferent starting values have been chosen. Only points for 
large enough n are included in the chain analysis such that 
the projections of each parameter distribution are consis- 
tent within statistics. 

While this is possible for cases with well-measured and 
thus strongly constraining obscrvables and a small num- 
ber of parameters (e. g. fits of a high-scale SUSY model 
like mSUGRA using LE and LHC measurements), this 



approach fails for more challenging problems like fits of 
an 18-dimensional more general MSSM. There are around 
20 million points in a combination of several Markov 
Chains with different starting points. These can not be 
checked efficiently for the effect shown in Figure 1 . Neither 
is it technically possible to provide sufficiently long chains 
due to computing limitations. Therefore, the Bayesian in- 
terpretation is used in this paper only in the case of fits 
of the mSUGRA model to the LHC data for illustration. 

The Frequentist interpretation of the Markov Chain is 
used as a default in this paper. It does not make use of 
the sampling density directly, but employs only the ob- 
tained x 2 values found in the chain for each parameter 
point. The best fit point is directly defined by the pa- 
rameter point with the lowest \ 2 = Xmin ( or equivalently 
the point with the largest likelihood £ = £ max )- To ob- 
tain uncertainties for a subset of parameters this approach 
scans over all parameters except for those under study 
and chooses the scanned parameters such that C is max- 
imised for each point in the studied parameter subspace. 
This procedure yields a profile likelihood. The la {2 a) 
uncertainties in the one-dimensional case arc defined by 
Ax 2 = X 2 — Xmin = 1 (4) an( l, in the two-dimensional 
case, by Ax 2 = X - Xmm = 2.3(5.99) [86]. In the limit 
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Fig. 2: Frequentist interpretation of Markov chain fitting 
mSUGRA using only LE measurements for two different 
starting points. 



of Gaussian parameter distributions, i. e. a locally linear 
relation between measurements and parameters, the 1 a 
environment for a given set of measurements covers an 
area which contains the true parameter point in 68 % of 
all possible experimental outcomes. 

This approach has several advantages over the 
Baycsian interpretation. First, it does not depend on prior 
distributions, since the likelihood in the hidden dimen- 
sions is not integrated. Instead the Markov Chain is sim- 
ply used an efficient scanning technique for the parameter 
space: for each bin in the histogram of the paramcter(s) 
under study the point with the lowest x 2 is chosen. Sec- 
ond, it is sufficient to scan different regions in parameter 
space with different granularity, as long as the obtained 
point density is fine enough to derive a smooth surface. 
This approach requires significantly less points than the 
Bayesian approach. 

In order to test if the sampling density is sufficient 
to derive a smooth x 2 surface, the following procedure 
is used. For each result presented in this paper, at least 
two different parameter points in the parameter space are 
chosen as starting point. These starting points are chosen 
individually for each fit after an initial Markov Chain run 
such that they lie approximately 2 a away from the esti- 
mated best fit point. Then, Markov Chains are computed 
for each of the two starting points. An example for the 
comparison of the two results can be found in Figure 2. 
For a fit to be accepted, it is required that the differences 
in the shape and area of the 1 a and 2 a regions are within 
the differences of the binning, which is derived for each 
chain separately by dividing the observed 2 a region into 
25 x 25 bins. This procedure to ensure the robustness of 
the results turned out to be very important, in particular 
in the case of more complex fits, e. g. the MSSM fit to LE 
and LHC measurements and the mSUGRA fit to LE mea- 
surements. In these cases, several million sampling points 
turned out to be required. 



The results of the Bayesian and Frequentist interpreta- 
tions are compared for the case of mSUGRA fitted to LHC 
measurements. For more complex fits there is not enough 
confidence in the results of Bayesian interpretation for the 
reasons explained above (see also [85]). 

As described above, the uncertainties on the model 
parameters can be derived from the A\ 2 values of each 
parameter point. It is important to note that the identifi- 
cation of Ax 2 = 1 with the 68 % uncertainty region is tech- 
nically only true for problems with Gaussian uncertainties 
of the observables and linear dependencies between pa- 
rameters and observables, hence for Gaussian parameter 
distributions. As visible already from Figure 2, this as- 
sumption is not fulfilled for all fits. Furthermore, for the 
interpretation of the fit result and especially for the deriva- 
tion of conclusions on RGE running of parameters or the 
prediction of other observables which are not yet used in 
the fit, it is important to also derive correlations between 
the model parameters. The derivation of linear correla- 
tion coefficients can be achieved by interpreting the ob- 
tained x" 2 surface of each combination of two parameters 
of the model as a numerical version of the Hesse-matrix of 
the problem. Normally, e. g. in MINUIT, the Hesse ma- 
trix is of course an approximation of the true x 2 surface. 
Hence the correlation coefficient can be calculated from 
the two-dimensional histogram of —Ax 2 for each param- 
eter combination. For practical reasons the histogram is 
constrained to the region of Ax 2 < 5.99, corresponding to 
the two-dimensional 95 % uncertainty region around the 
best fit point. This is justified because the linear correla- 
tion coefficient are an approximation of the full non-linear 
correlations and are dominated by the area around the 
minimum. 

The Frequentist analysis of the parameter space using 
Markov Chains could be further refined using MINUIT 
around the x 2 minimum found in the Markov Chains in 
order to better determine the position of the absolute min- 
imum, independent of binning effects. However, in prac- 
tice each set of Markov Chains with sufficient sampling of 
the Ax 2 < 5.99 region, as required above, ensures a suffi- 
ciently fine sampling around the absolute minimum, such 
that no significant improvement of a further refinement of 
the minimum in a MINUIT minimisation can be seen. 

While Markov Chains provide a powerful tool for the 
study of complex parameter fits, there are a few short- 
comings which make it desirable to cross-check the results 
with an alternative technique and to provide solutions for 
tasks which cannot be solved by Markov Chain Monte 
Carlos. First, it would be desirable to use a method which 
is independent of both the problem of priors (as the case 
in the Bayesian interpretation) or assumptions on (local) 
linearity of measurements and parameters (as the case in 
the Frequentist interpretation). 

Second, the Markov Chains cannot be used to study 
how well a certain model can be distinguished from an- 
other model by the data, or how one interpretation of the 
observables is distinguished from another interpretation 
within the same model. Markov Chains just supply the 
P-value of the best fit point of each model or each inter- 
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prctation. They do not provide the probability to get a 
better P-value for the "wrong" model than for the "cor- 
rect" model under the assumption that one of the models 
is realised. 



3.1.2 Toy Fits and Simulated Annealing 

In addition to the Markov Chain analysis, Toy Fits are 
used to obtain an independent estimate of the parame- 
ter uncertainties and to compare different models or data 
interpretations quantitatively. The Toy Fit analysis con- 
sists of two steps. First, Markov Chains or minimisation 
through Simulated Annealing (for a description of the im- 
plemented algorithm of Simulated Annealing see [17]) is 
used to find the absolute x 2 minimum of a given problem. 
Second, Monte Carlo Toy data are created around the ob- 
servables corresponding to best fit point. In the first step, 
the best fit parameters P are determined. Then the set 
of observables M a corresponding to this set is calculated. 
This set is then used to create N different MC Toy sets of 
pseudo-measurements Mi (i. c. other possible experimen- 
tal outcomes which are consistent with that parameter 
set) by smearing around M Q according to the Gaussian 
uncertainties and correlations as defined in covm- For each 
of the N MC Toy sets M", a fit is performed using Sim- 
ulated Annealing followed by a MINUIT minimisation at 
the minimum of the Simulated Annealing fit. This proce- 
dure yields a set of N "best fit" parameter points Pi . The 
distributions and correlations of Pi are then interpreted 
as the expected distributions of all possible experimental 
outcomes given the best fit parameter set P . 

In contrast to the Frequentist interpretation of the 
Markov Chain Monte Carlo, the uncertainties and corre- 
lations can be directly calculated from the (co-)variances 
of Pi. Therefore, the results represent an estimate for the 
expected distribution of all possible results (including all 
possible non-Gaussianities). They even include possible 
secondary x 2 minima, which are turned into the absolute 
minimum for a subset of the observable set Mi , where sta- 
tistical and systematic uncertainties of the measurements 
can invert the order of the different x 2 minima. 

Toy Fits allow for a robust cross-check of the validity 
of the results: For problems which arc sufficiently Gaus- 
sian the distribution of the Xmin va hies has to be consis- 
tent with a x 2 distribution for n — m degrees of freedom, 
where n is the number of observables and m is the num- 
ber of parameters. This criterion is checked by fitting the 
X 2 distribution to each Xmin histogram and requiring that 
the fitted number of degrees of freedom agrees with the 
expectation for the problem within 2 a for a fit to be ac- 
cepted. 

In addition, as for Markov Chains, starting values for 
the parameters have been varied and the result is required 
to be independent of the starting value up to statistical 
fluctuations. 

However, there is also a disadvantage of the Toy Fits 
with respect to the Markov Chains. In case of two different 
X 2 minima very far apart from each other (with respect 
to the parameter uncertainties, and with a very high x 2 



barrier in between) which both have almost identical min- 
imal x 2 1 one se t °f Toy Fits with one given starting point 
_P s i of the fits would not necessarily find both minima, re- 
quiring possibly more than two different starting points. 
In the limit of infinite computing power this can be over- 
come by starting the individual Toy Fits from randomised 
positions within the parameter space in order to scan for 
additional minima. If additional minima with acceptable 
P-values are found and if the x 2 barriers between them are 
too high to be traversed within one minimisation, different 
sets of Toy Fits with starting values around the different 
minima can be directly treated as different models using 
the prescription proposed in Section 3.2. This also allows 
to assign a numerical value to the parameter uncertainties 
stemming from the different minima. 



3.2 Model Discrimination and Ambiguities 

All techniques discussed above are directed towards deter- 
mining the uncertainties of the parameters of one model in 
the presence of a given set of data with one fixed interpre- 
tation of the data. In this paper a new method is proposed 
to use Toy Fits to solve the two remaining tasks, namely to 
measure how often fluctuations of the data lead to wrongly 
identifying a different model than the true model as the 
model with the best P-value, and determining how ambi- 
guities in the mapping between measurements (i.e. edges 
in LHC mass spectra) and predicted observables, translate 
into the uncertainties of the model parameters. 

The first of these problems cannot be answered di- 
rectly by using the 'P-valuc. It represents the probability 
that the given experimental data are observed given a cer- 
tain model with certain best fit parameters, and hence can 
be used to determine which models (or which interpreta- 
tions of the data in a given model) are acceptable and 
which models are rejected by the data due to their low 
P-valuc. However, it cannot be used to determine how 
likely it is to obtain a certain P-value for a given model 
if the data are actually caused by another model. Such 
a situation can occur often for not very well constrained 
problems, as e. g. visible in Sections 4.1.1, 4.1.3 and 4.3.1. 
This question can be answered by the Toy Fits directly by 
fitting an arbitrary number of different models (or differ- 
ent interpretations of the data in the same model) j with 
different predictions 0{P\) to the same set of smeared 
measurements Mi. Then the probability p w to prefer the 
"wrong" models over the "correct" model, from which the 
observable set M is derived, can be directly determined 
by counting how often one of the "wrong" models achieved 
a better Xmin than the "correct" model. Since in reality 
the "correct" model is not known, this interpretation can 
be repeated for each model which yields a reasonable V- 
valuc. Examples using this procedure are presented in Sec- 
tion 4.3.1. 

The second problem can be solved directly using the 
simultaneous Toy Fit method proposed above. The uncer- 
tainty of a model parameter is defined and measured as the 
square root of the variance of the parameter distribution. 
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Each entry Pi in the Toy Fit parameter distributions cor- 
responds to one individual best fit point for one model, 
one interpretation of the data and one individual set of 
smeared measurements Mi. 

In the presence of several interpretations of the data 
(or e. g. different possible values of a discrete parameter 
of the model, i. c. in all cases where different model/data 
combination can be fitted with the same parameters and 
the same measurements) the distributions of the Pi can 
be exchanged against the distribution of the Pj, where 
for each smeared simulated Toy measurement i the model 
or interpretation j yielding the best fit is chosen which 
yields the best x 2 . This is the natural extension of the 
method exploited usually in Toy Fits and described in 
Section 3.1.2, since in the presence of one given set of 
measurements Mi all remaining ambiguities of the inter- 
pretation of the data would be tested and for the final 
result the interpretation with the lowest x 2 (i- e - highest 
P- value) would be chosen. The uncertainty stemming from 
the ambiguities is then taken into account by not citing 
the parameter uncertainties of the best fit interpretation of 
the data only, but by citing the uncertainty including the 
possibility that a different interpretation would have been 
chosen as the one with the best P-value, as proposed here. 
Note that this interpretation naturally leaves the parame- 
ter distribution Pi unchanged if only one model is always 
yielding the best 'P-valuc, i. e. if the method proposed 
above determines that there is % probability to prefer a 
"wrong" model over the "correct" model. An application 
of this method is presented in Section 4.3.1. 



4 Results 

In this chapter, results are presented in terms of allowed 
areas in the SUSY parameter space. Also, allowed regions 
for SUSY particle masses are calculated from the fitted 
parameters. 

In Section 4.1, the high-scale SUSY models mSUGRA 
and GMSB are tested against presently available mea- 
surements. Predictions for discovery at the LHC and the 
expected range of SUSY masses are presented. In Sec- 
tion 4.2, the constraint from LHC observables are stud- 
ied alone. In Section 4.3, the expected measurements at 
the LHC are combined with already available observables. 
Here both mSUGRA and a more general 18-parameter 
MSSM are studied. Finally, in Section 4.4, the impact of 
SUSY precision measurements at the ILC together with 
£ mt = 300 fb -1 of LHC data is studied, in particular for 
the general MSSM. 



4.1 Present Low-Energy Observables 

Several studies of the mSUGRA parameter space in the 
light of different sets of available measurements have been 
performed recently [44,85]. In this section, parts of these 
studies are repeated and extended. Emphasis is laid on 
understanding the impact of the most important observ- 
ables on the parameter uncertainties, and on the study 



Table 3: Result of the fit of the mSUGRA model with 
sign(/x) = +1 including four additional SM parameters 
to all measurements listed in Table 1. The minimum % 2 
value is 20.6 for 22 degrees of freedom, corresponding to 
a-P-value of 54.4%. 



Parameter 


Best Fit 




Uncertainty 


sign(/j) 


+ 1 






a s 


0.1177 


± 


0.0020 




127.924 


± 


0.017 


m z (GeV) 


91.1871 


± 


0.0020 


mt (GeV) 


172.4 


± 


1.09 


tan f3 


13.2 


± 


7.2 


M 1/2 (GeV) 


331.5 


± 


86.6 


Mo (GeV) 


76.2 




+79.2 
-29.1 


A (GeV) 


383.8 


± 


647 



of the uncertainties of SUSY particle mass predictions. 
The methods outlined in Section 3 are used to extract pa- 
rameter correlations. The impact of the SM parameters is 
also studied. In addition, a comparison of the predicted 
mSUGRA sparticle spectrum with a GMSB spectrum is 
presented. 

In order to study the effect of the different observ- 
ables, first a baseline fit with all observables from Table 1 
is performed. In order to accommodate the uncertainties 
of the most important observables, (<? — 2) M in terms of its 
SM prediction, and ficDMh 2 in terms of its origin, several 
other fits are then shown. Additionally, these observables 
are treated in different ways, removed one by one or re- 
moved simultaneously 



4.1.1 Fits of the mSUGRA Model 

The mSUGRA (or sometimes the CMSSM, which has n 
as additional free continuous SUSY parameter) scenario 
is the best studied SUSY scenario at colliders. For the 
purpose of this paper, it is appealing to study it in detail 
using the low energy observables, since experimental stud- 
ies of possible LHC measurements are available from AT- 
LAS [87] and CMS [88] for several well studied mSUGRA 
parameter points. Knowing the allowed parameter space 
using the already existing observables allows the selection 
of a mSUGRA scenario studied at LHC which is within 
the currently experimentally allowed parameter region. 

The baseline fit is performed using all measurements 
from Table 1, i. e. assuming that the current predictions 
for the SM contribution to (g — 2) M arc correct and that 
the cold dark matter in the Universe is entirely caused by 
the relic density of SUSY LSPs, which in the mSUGRA 
scenario has to be the lightest neutralino Xi- For all fits 
in this paper, if not stated otherwise, all parameter points 
with charged stable LSPs are excluded from the fit. The fit 
has been performed using the Markov Chain Monte Carlo 
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technique described in Section 3.1.1 and applying the Fre- 
qucntist interpretation. Using two different starting values 
of the Markov Chains, approximately 20 million points 
per chain have been tested with consistent results for the 
different starting points. The sampling density around the 
minimum is large enough such that an additional MINUIT 
minimisation is not necessary. The x 2 minimum of 20.6 
for 22 degrees of freedom, corresponding to a P-value of 
54.4%, can be determined with good precision from the 
distribution of the sampled points. The P-valuc is signifi- 
cantly larger than the P-value of the SM fit [89] of around 
18 %, which does not include the limits from direct Higgs 
searches. With the direct Higgs search limit included, the 
P-value is 17% [45]. However the analysis of the variable 
pulls explained below shows that this is not directly due to 
a better description of the SM precision variables, which 
cause the moderate to low P- value of the SM fit. It can 
be seen that mSUGRA provides an excellent description 
of the presently available precision data. 

Correlations among the obscrvablcs of Table 1 , as de- 
scribed in Section 2.1.3, are studied in a separate fit. While 
the minimal x 2 is lowered by 1.1 due to the correlations, 
the results of the allowed parameter regions for the fit in- 
cluding correlations are identical to those of the baseline 
fit without correlations in the input observables. The rea- 
son for that lies in the fact that (g — 2) M and ^cdm^ 2 
(which are uncorrelated) constrain the parameter space 
more strongly than any other combination of variables 
(see Section 4.1.2). Therefore there is no significant im- 
pact of the correlations among the electroweak precision 
observables. 

The result of the baseline fit is given in Table 3. Due 
to the deviation of (g — 2)^ from the SM prediction to- 
wards larger values, a positive sign of fi is preferred, which 
is chosen for this fit. A comparison with sign(/i) = — 1 
is shown below. It can be seen that moderate values of 
tan/3 between 5 and 20 are preferred. The gaugino mass 
parameter M 1 / 2 is expected in the range of 200 GeV to 
400 GeV, while a low scalar mass parameter Mo between 
50 GeV and 150 GeV is preferred. The universal trilinear 
coupling Aq is not very well constrained and is expected 
to be between -300 GeV and 1000 GeV. As it is shown 
below, this region of parameter space is very favourable 
for early discoveries at LHC and for a rich phenomenol- 
ogy at ILC. In Table 3, as well as in the following tables, 
we quote symmetrical uncertainties whenever upper and 
lower uncertainties agree within 20%, and asymmetrical 
uncertainties are given otherwise. It should be noted that 
the probability densities for the fitted parameters are usu- 
ally not Gaussian, i. e. that the one-dimensional 2 a un- 
certainty at Ax 2 = 4 are not twice the 1 a uncertainty at 
A X 2 = 1. 

Note that the uncertainty of the SM parameters ex- 
actly corresponds to the uncertainties of the measure- 
ments to which the given parameter is 100% correlated 
(see Table 1). This means that the direct measurement of 
the parameters is so precise that the additional observ- 
ables from Table 1 do not play any role in their determi- 



nation, which has consequences for their correlations with 
other parameters. 

It also has to be noted that Gp is omitted from Ta- 
ble 3, since Gf is measured so precisely that virtually no 
effect of the inclusion of Gf into the fit can be observed. 
This has been checked using a second baseline fit including 
Gf, which yields identical results for the other SM param- 
eters and the mSUGRA parameters. Hence Gf has been 
omitted from all subsequent fits to save computing effort. 
For completeness however, it is included in the discussion 
of the parameter correlations. 

In addition to the parameter uncertainties, it is im- 
portant to study the parameter correlations. Calculating 
uncertainties of quantities derived from the parameters 
docs not only depend on the absolute value of the un- 
certainties, but also on the correlations. The prescription 
proposed for the calculation of correlations from Markov 
Chain Monte Carlos in Section 3.1.1 is used for the re- 
sults in Table 4. It is interesting to note that while there 
are significant correlations among the mSUGRA parame- 
ters, there is no correlation in excess of 10 % between any 
SM parameter and any mSUGRA parameter. In addition, 
there are hardly any correlations among the SM param- 
eters themselves. The reason for this is the fact that ev- 
ery SM parameter can be measured using an observable 
to which the parameter is correlated 100%. This ensures 
that each SM parameter is fixed by itself to a strong pre- 
cision (due to the impressive success of the SM precision 
measurements) without any impact of any other param- 
eter. It should be noted however, that the correlations 
of Table 4 are the linear correlations of the parameters 
within the two-dimensional 2 a area. For some choices of 
observables (e. g. without the use of (g — 2)^, see below) 
the preferred parameter space exhibits strong non-linear 
correlations, restricting the reliability of this method to 
well-constrained fits like the baseline fit described here. 

The individual pull of each observable is shown in Fig- 
ure 3. Again it can be seen that the SM parameters can 
all be fitted individually exactly to the values of their 
corresponding observables. When comparing the result 
with the SM fits in [89,45], the main difference is that 
m/j is not a free parameter anymore, but a function of 
the SUSY parameters. This has the effect that the elec- 
troweak precision observables can be satisfied with a sig- 
nificantly larger P-value than for the SM fit. The pre- 
ferred value of nih is 113.3 GeV and this is below the 
95 % CL lower limit of the SM Higgs boson searches, which 
are applicable to mSUGRA, as outlined in Section 2.1.4. 
The preferred value is within 1/3 of 1 a of the theoreti- 
cal uncertainty with respect to the 95% CL lower limit 
at m/j > 114.4 GeV, and even within approximately 1 a 
of the lowest mass point where the CL S of the combined 
searches is 0.5 at around rrih = 116.5 GeV. The latter is 
corresponding to the mean value of a "measurement" of 
77i/i, if the searches are interpreted as such [54]. 

It can be observed that apart from the change in rrih, 
the fit results are remarkably similar to the SM fit. This 
is due to the decoupling of the sparticles heavier than the 
electro-weak scale from the processes contributing to the 
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Table 4: Correlations of the mSUGRA model and the SM parameters (with sign(^) = +1) using all measurements 
listed in Table 1. Very small correlations between SM and mSUGRA parameters is observed while there is significant 
correlation among the mSUGRA parameters. 



Parameter 


a s Gf 


1/aem 


mz 


mt 


tan f3 


Mi/ 2 


Mo 


Ao 


Ois 


1.000 -0.005 


0.006 


-0.003 


-0.007 


-0.003 


-0.005 


-0.013 


-0.009 


Gf 


1.000 


-0.003 


-0.001 


-0.003 


-0.002 


-0.022 


0.007 


-0.006 


1/otem. 




1.000 


-0.008 


0.006 


-0.001 


-0.003 


0.009 


0.007 


m z 






1.000 


-0.0348 


0.053 


0.035 


0.046 


0.029 


m t 








1.000 


0.075 


0.088 


0.075 


0.093 


tan j3 










1.000 


0.358 


0.833 


0.457 


Mi /2 












1.000 


0.449 


0.236 


M 














1.000 


0.632 


A 
















1.000 



mSUGRA fit to LE 



m t 


172.4 ± 1.2 


172.4 


m b 


4.2 ±0.1 7 


4.2 


m z 


91.1875 ±0.0021 


91.1871 


a. 


0.1 176 ±0.0020 


0.1177 


G F 


1.16637 10" 5 ±10' 1 ° 


1.1663710 


"em 


127.925 ±0.016 


127.924 


m h > 


114.4 


113.3 




41 .54 ±0.04 


41.48 


ir 


0.01714 ±0.00095 


0.01644 


A t 


0.1 465 ±0.0032 


0.1480 


A, 


0.1 51 3 ±0.0021 


0.1480 


A c 


0.67 ± 0.027 


0.67 




0.923 ± 0.02 


0.935 




0.0707 ± 0.0035 


0.0742 


0.0992 ±0.0016 


0.1038 


Rc 


0.1721 ±0.003 


0.1722 


R b 


0.21 629 ±0.00066 


0.21604 


R i 


20.767 ± 0.025 


20.746 


r z 


2495.2 ± 2.51 


2495.1 


sine sff 


0.2324 ±0.0012 


0.2314 




80.399 ± 0.027 


80.380 


"dm 
(9-2), 
BR(b— > sy) 


0.1 099 ±0.01 35 


0.1115 


3.02 10" 9 ± 9.0 10' 1 ° 


2.55 10' 9 


1.117± 0.122 


1.009 


BR(b-> xv) 


1.15 ±0.4 


0.96 


BR(B s -> X s ll) 


0.99 ± 0.32 


0.99 


BR(K% Iv) 


1.008 ±0.014 


1.000 




0.92 ±0.14 


1.03 


A(m s ) 


1.11 ±0.32 


1.03 


Am/ A m 


1.09 ±0.16 


1.00 




12 3 
(Meas.-Fit)/ o 

Fig. 3: Pull for low energy observables used in the 
mSUGRA parameter fit using the observables from Ta- 
ble 1 and the best fit point from Table 3. 



precision measurements, hence rendering the result very 
SM-like. It can be seen that while rendering rrih naturally 
heavier, mSUGRA does not mitigate the tension in some 
of the SM precision observables. 

Apart from looking at the numerical results of the 
baseline mSUGRA fit, it is interesting to study the allowed 
parameter regions in detailed two-dimensional projections 
of each mSUGRA parameter against another parameter, 



as shown in the profile likelihood plots obtained using the 
Frcquentist interpretation of the Markov Chains in Fig- 
ure 4. The SM parameters have been suppressed due to 
their negligible correlations with themselves and with the 
mSUGRA parameters (see Table 4). The plots show the 
two-dimensional 95 % CL allowed region of the fits in blue 
(corresponding to A\ 2 — 5.99) and the one-dimensional 
68% region in red (corresponding to A\ 2 = 1). There 
are two reasons why the unusual choice of showing the 
A\ 2 = 1 curve in a two-dimensional plot has been made. 
First because by this the one-dimensional uncertainties 
of the parameters can be directly read off from the plot, 
which is not possible for the choice of the two-dimensional 
68% CL area at A\ 2 = 2.3. Second, the two-dimensional 
95 % CL area gives a good indication of the experimentally 
allowed area, while two-dimensional 68 % CL area leaves 
a large room for parameter points outside the 68 % con- 
tour. The most common projection is shown in the upper 
left plot, which compares the allowed region in Mq and 
Mj/2 between the baseline fit described above and a fit 
requiring only il t heoh 2 < f2 b s h 2 , as described below in 
Section 4.1.2. It can be seen that on the upper left side 
the allowed region is directly adjacent to the excluded re- 
gion where is the LSP, which is excluded from the fit 
because a stable charged LSP is in conflict with cosmo- 
logical measurements. The results also show that on the 
2 cr level M 1/2 < 800 GeV and M < 600 GeV is expected, 
yielding relatively light sparticles and hence good discov- 
ery prospects for the LHC. The upper right plot shows 
the projection of tan/3 versus Mo, showing that on the 
2 a level tan /3 values are allowed between 3 and 50. This 
means that in a full simultaneous fit of all parameters 
tan P is not constrained strongly. This implies that a wide 
variety of different signatures is possible, as outlined in 
more detail below. In the middle plot on the right side it 
can be seen again that Aq is not constrained very strongly, 
leaving a lot of room for the phenomenology of the third 
generation. However, this has only a small effect on the 
discovery prospects at LHC, since both tan/3 and A have 
limited effects on the first two generations and the gaugi- 
nos, which are most important for the discovery modes at 
the LHC. 
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Fig. 4: mSUGRA parameter regions compatible with relic density constraint and with all low energy measurements 
for various parameter combinations and sign(/i) fixed to +1. For both cases two-dimensional 95 % confidence level 
and one-dimensional 68 % confidence regions arc shown. From the latter 1 a uncertainties for individual parameters 
can be derived from a projection of the area to the respective axis. The dashed lines represent the result when all 
observables are included. For the full lines, the measured relic dark matter density is only regarded as an upper limit 
of the predicted SUSY dark matter density. All other observables remain the same. 
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Fig. 6: SUSY mass spectrum as predicted by mSUGRA 
parameter fit to low energy measurements with sign(/i) 
fixed to +1. 



Fig. 5: Allowed mSUGRA parameter space from Table 3 
overlaid upon the expected ATLAS discovery reach using 
1 ftT 1 of data at ,/s = 14 TeV [87]. 

The direct comparison of the allowed parameter space 
of the baseline fit with the expected discovery reach of 
the ATLAS experiment is shown in Figure 5. The AT- 
LAS discovery reach plot is calculated for tan (3 = 10 and 
Aq = GeV. Although the central values of the baseline 
fit of these parameters do not exactly agree with these 
settings, it is justified to compare the fit results with the 
ATLAS discovery reach plot. This is the case because first 
the fixed values in the discovery reach plot are within the 
uncertainties of the fit and second because the discovery 
reach depends mainly on the gluino and first generation 
squark masses. These are very insensitive to Aq and not 
very sensitive to tan (3 due to the absence of mixing on tree 
level. The lines in Figure 5 correspond to the boundaries 
of the 5 a discovery region. The most sensitive search is 
expected to be an inclusive measurement of the effective 
mass spectrum of 4 jets and missing transverse energy. 
It can be seen that almost the entire mSUGRA parame- 
ter space allowed at 95 % CL is observable already with 
1 fb _1 of well-understood data at = 14 TeV of the 
ATLAS experiment [87]. 

From the Markov Chain results, also the probabil- 
ity densities for the sparticle masses can be deduced. In 
Figure 6 the expected masses of the Higgs bosons and 
sparticles are shown for the baseline fit. The red lines 
indicate the masses corresponding to the best fit pa- 
rameter point. The dark (light) blue regions denote the 
one-dimensional 68 % (95 %) CL area, corresponding to 
the A\ 2 < 1 {A\ 2 < 4) region around the absolute 
minimum. The expected spectrum shows several distinct 
features: The Higgs boson mass is well constrained to 
= 113.3 ± 2.5 GeV, just above the LEP exclusion 
bound. The lightest SUSY particle (LSP) is the stable 
Xi at m^o = 130 ± 35 GeV with a relatively small mass 
difference to the next-to-lightest SUSY particle (NLSP) 



fi at mfj — 140 ± 25 GeV, with the mass difference be- 
tween NLSP and LSP being rather precisely constrained: 
77if 1 — m^a = 9.5 ± 2.5 GeV. If mSUGRA is realised in 
Nature, this would provide challenging experimental con- 
ditions for precision measurements at the LHC, since the 
very small mass difference between NLSP and LSP would 
provide many final states with very soft particles, domi- 
nated by r-leptons, as the last reconstructable particles in 
the SUSY decay chain. In comparison with the discovery 
reach from Figure 5 it can be seen that mSUGRA would 
allow for relatively easy discovery in inclusive channels 
but difficult precision measurements in exclusive recon- 
struction of SUSY cascades. 

Furthermore, all light neutralinos, charginos and slep- 
tons are expected well below m « 600 GeV. The squarks 
and the gluino could be found below 1.6 TeV. The heavy 
Higgs bosons are expected below m H /A < 1.2 TeV and 
would be challenging to discover even with more than 
lOftV 1 of luminosity at the LHC, since tan/3 is not large 
enough [87]. 

The baseline fit uses a fixed value of the discrete pa- 
rameter sign(/i) = +1, since this choice is preferred by 
the positive deviation of the measured anomalous mag- 
netic moment of the muon a° xp from the SM prediction. 
However, the difference in the overall agreement between 
the choice of sign(/j) = +1 and sign(/i) = — 1 has to be 
assessed with a global fit to the same observablcs as used 
by the baseline fit, but for sign(/i) = — 1. The result of the 
this fit (again using Markov Chains) is shown in Table 5. 
Due to the tension in a° xp — a„ of 3.4 cr, the minimal 
X 2 of the fit rises to 31.09, corresponding to a 'P-valuc 
of 9.4%, compared to 54.4% for the baseline fit. This re- 
sult shows that sign( J tt) = — 1 is disfavoured but not ex- 
cluded. The preferred parameter regions for both choices 
of sign(/i) = — 1 show similar preferred values for tan/3 
and M 1 / 2 - For sign(/i) = —1, much larger values of Mq 
are preferred. 
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Fig. 7: mSUGRA parameter regions compatible with the relic density constraint and with all low energy measurements 
for various parameter combinations and sign(/i) fixed to —1. For both cases two-dimensional 95 % confidence level 
and one-dimensional 68 % confidence regions are shown. From the latter la uncertainties for individual parameters 
can be derived from a projection of the area to the respective axis. 



Table 5: Result of the fit of the mSUGRA model with 
sign(/i) = — 1 including four additional SM parameters 
to all measurements listed in Table 1. The minimum % 2 
value is 31.1 for 22 degrees of freedom, corresponding to 
a^-value of 9.4%. 



Parameter 


Best Fit 




Uncertainty 


sign(^) 


-1 






a s 


0.1177 


± 


0.0018 


l/a em 


127.924 


± 


0.014 


m z (GeV) 


91.188 


± 


0.0022 


m t (GeV) 


172.5 


± 


1.02 


tan f3 


9.6 




+ 17.8 

-4.5 


M 1/2 (GeV) 


125.1 




+70.0 
-25.1 


M (GeV) 


2313.5 




+622 
-940 


A (GeV) 


-29.1 


± 


2048 
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Fig. 8: SUSY mass spectrum as predicted by mSUGRA 
parameter fit to low energy measurements with sign(ju) 
fixed to —1. 



The two-dimensional projections of the mSUGRA pa- 
rameters for sign(/i) = —1 are given in Figure 7. In 
contrast to the baseline fit, where one continuous two- 
dimensional 95 % CL area is observed, two distinct areas 
are found for sign(^i) = —1. One is located at large -Mi/ 2 
and small Mo, while the other one exhibits the opposite 
signature. The reason is that in both regions the negative 
SUSY contribution to (<? — 2)^ is reduced, while a correct 
contribution to the cold dark matter relic density is pre- 
served. The region of both large M 1 / 2 and Mo is therefore 
cut out by the constraint of I2cdm^ 2 - 



As before, the result can also be expressed in terms of 
the preferred sparticle masses. Fig 8 shows the best fit and 
the 68 % and 95 % CL areas of the observable masses for 
the mSUGRA fit with sign(^) = — 1. Since large Mo and 
small My j 2 is preferred over large M^ and small Mo, as 
visible in Figure 7, the expected gaugino masses in this 
scenario are small. With m^o = 55^10 GeV a very light 
LSP is predicted just above the LEP mSUGRA limit of 
m^o > 47 GeV [90]. In contrast to the baseline fit, heavy 
sleptons and squarks are expected due to the large value of 
Mq. The gluino is predicted to be lighter than all squarks 
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and sleptons such that it decays via three-body decays to 
two jets and the LSP. 

No tension between the predicted lightest Higgs boson 
mass with the LEP limits is present for sign(/i) = — 1. 
The obtained Higgs mass value is rrih = 118.3 ± 3GeV 
well above the LEP bound at rrih > 114.4, GeV and also 
above the preferred LEP Higgs mass of rrih ~ 116.5 GeV. 

In addition to the mass spectrum and the important 
mass differences, the expected dominant sparticle decay 
modes can be studied. A selection of branching fractions 
of the mSUGRA parameter points preferred by the base- 
line fit is shown in Table 6. The central values correspond 
to the best fit, the uncertainties are given by A\ 2 < 1 
with respect to the best fit. It can be generally observed 
that the uncertainties of the expected branching fractions 
are very large, often close to 100%. Also, there are only 
few sparticles with only one relevant decay mode. There- 
fore a rich phenomenology with many competing decay 
modes can be expected at the LHC. Branching fractions 
of the electro-weak gauginos are typically largest for the 
third generation, with smaller contributions from the first 
and second generation. In connection with the small mass 
difference between and to^o , the larger branching frac- 
tions into r leptons compared to other leptons represent 
a challenge for the measurement of ratios of branching 
fractions and di-tau mass endpoints. 

Following the expected decay chains, the gluino decays 
into squarks with a slight but not dominant preference for 
bi/2- Generally, the branching fractions into the spartners 
of the right-handed degrees of freedom are larger by al- 
most a factor of 2 than the branching fractions into the 
left-handed counterparts. The right squark q~R decays into 
Xi<7 almost exclusively. As expected, the decay of the q~L 
is predicted to be more complex with decays into xW and 
Xi q' ■ The chargino has a preference for decays into t\v t , 
but includes significant contributions of decays into sneu- 
trinos and into y§W • The x.2 nas a similarly rich spec- 
trum of decays, with a dominance of fir, but with small 
contributions of sneutrinos, other leptons and Xi^° an( i 
XiZ°. The sleptons exhibit branching fractions of close to 
100 % into X? and the corresponding lepton. 

The results of the baseline mSUGRA fits with 
sign(/i) = +1 and —1 provide clear predictions for the 
expected measurements at the LHC. However, there are 
several interesting questions which remain. First, it re- 
mains to be assessed which measured observables con- 
strain the parameter space and the regions of LHC ob- 
servables in which way, i. e. which features of the mea- 
surements dominate the prediction of parameters and fu- 
ture observables. Second, as outlined in Section 2.1, the 
interpretation of some of the measurements in terms of 
mSUGRA is not necessarily unique. Third, alternative 
SUSY breaking models may predict other features for the 
LHC. These questions are addressed in the next sections. 



Table 6: Expected branching fractions of the SUSY par- 
ticles in the mSUGRA model with sign(/i) = +1. The 
results shown for the first generation are also valid for the 
second generation. 



Decay Mode Expected Branching Fraction Uncertainty 



X°2 -> TlT 


0.46 


+0.38 
-0.44 


X.2 &T! V r 


0.076 


+0.039 
-0.067 


X% -*■ e R e 


0.040 


+0.044 
-0.038 


X°2 - X?^ 


0.036 


+0.13 

-0.035 


X°2 - x\z 


0.018 


+0.0098 
-0.018 


X2 -> e L e 


0.00018 


+0.14 
-0.00018 


X2 -*■ f 2 r 


0. 


+0.014 
-0 


Xt -» hV-r 


0.40 


+0.42 
-0.39 


xt -> Ve L e 


0.15 


+0.10 
-0.15 


xf -> &T1 t 


0.15 


+0.10 
-0.15 


xt - i&w* 


0.12 


+0.079 

-0.12 


Xt - > 





+0.14 

-0 


g — > ulu 


0.052 


+0.020 
-0.023 


g — > u R u 


0.094 


+0.030 
-0.054 


g — > did 


0.051 


+0.018 
-0.021 


g — > d R d 


0.093 


+0.031 
-0.055 


g — > tit 


0.090 


+0.166 
-0.066 


g -> t 2 t 





+0.056 
-0 


g bxb 


0.179 


+0.022 
-0.056 


g — > b 2 b 


0.11 


+0.009 
-0.038 


e-R -»■ X?e 


1 


+0 
-0 


eL -»■ X?e 


1 


+0 

-0.23 


e-L -^xle 





+0.085 
-0 


h - x'h 


1 


+0 
-0 


T2 -> X° T 


0.99 


+0.011 
-0.42 


T% -> X2 T 


0.0007 


+0.040 
-0.0007 


T2 -> Xt V -r 


0.0044 


+0.078 
-0.0043 


UL XlU 


0.008 


+0.0049 
-0.0026 


UL -»• X°U 


0.31 


+0.013 
-0.0061 


UL -»■ X.t d 


0.64 


+0.019 
-0.0094 


UL X2 d 


0.025 


+0.0092 
-0.018 


UR -> X°iU 


0.99 


+0.012 
-0.0080 


ur -> xiu 


0.0078 


+0.0064 
-0.0072 


h -> x°t 


0.18 


+0.073 
-0.037 


tl -» X2t 


0.15 


+0.006 
-0.069 


ti -» xfb 


0.51 


+0.21 
-0.28 


ti -» xtb 


0.15 


+0.11 
-0.15 


bi ^ xtt 


0.37 


+0.13 
-0.039 


bi j&b 


0.25 


+0.12 
-0.068 


61 -> hw ± 


0.043 


+0.29 
-0.043 
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4.1.2 Fits of mSUGRA to Low-Energy Observables with 
reduced Sets of Observables 

In this section, the baseline fit with sign(/x) = +1 is mod- 
ified by fitting the same parameters to reduced sets of 
observables. As described in Section 2.1, both i^cDM^ 2 
and (g — 2)^ suffer from uncertainties concerning their in- 
terpretation in terms of SUSY. In this section we show 
that these observables provide the strongest constraints 
in the parameter space, therefore the uncertainty in the 
interpretation of these observables has to be evaluated for 
the predicted parameter space and collider observables. 

For i?cDM^ 2 there is very little doubt about the mea- 
surement itself. However, the cosmological measurement 
of the cold dark matter relic density does not imply that 
the SUSY LSP is solely responsible for the dark matter. 
Therefore, we study three different possibility: First, a 
stable neutral SUSY LSP in the context of a i?-parity 
conserving model is the only source of ^?cdm^ 2 , and the 
process of the LSP production after the Big Bang and 
the freeze-out of the LSP is completely understood, as 
e. g. implemented in [79]. This is assumed for the baseline 
fit in Section 4.1.1. Second, the SUSY LSP contributes 
to dark matter, but other unknown additional sources of 
dark matter are not excluded. This scenario is tested in 
this section by requiring that the LSP is stable and neu- 
tral and that the predicted cold dark matter relic density is 
smaller or equal than the observed one i? t heo^ 2 < flobsh 2 . 
This still includes the assumptions that the mechanisms of 
CDM creation are understood. Third, additional features 
like a CDM creation different from the current under- 
standing, or a meta-stable LSP, could cause differences be- 
tween the measured CDM relic density and the predicted 
CDM relic density. The maximal effect of such uncertain- 
ties are tested by removing i?cDM^ 2 from the observables 
in the fit. 

Similarly, three different scenarios can be distinguished 
concerning (g — 2) M , where there is an ongoing debate 
about uncertainties on the SM prediction (see e. g. [91,92, 
93,94]). Therefore, in addition to the baseline fit, which 
uses the current mean value of measurements and the SM 
and SUSY prediction from [53], deducing the hadronic 
corrections from e + e~ collision data, three other options 
are tested. First, in order to show the importance of the 
(g — 2) M measurement, it is removed from the observables 
used in the fit. Second, it is assumed for illustration that 
the current deviation between SM prediction and measure- 
ment is a statistical deviation or a insufficiency in the SM 
prediction, i. e. that there is no visible SUSY contribution 
to (g — 2) p. Third, the prediction deducing the hadronic 
corrections from r decay data [94] is used. Several com- 
binations of the above mentioned possibilities are studied 
in the following. 

Since the composition of the cosmos and its dynamics 
are a very dynamic field of study, it is a strong assump- 
tion that only the LSP contributes to cold dark matter. 
Therefore, a fit of mSUGRA with the preferred choice of 
sign(/x) = +1 is performed requiring J7 t heo^ 2 < ^obsh 2 
and leaving all other observables unchanged. In this fit, 
^theoft 2 docs not contribute to the \ 2 f° r ^thco^ 2 < 



Table 7: Result of the fit of the mSUGRA model with 
sign(/i) = +1 including four additional SM parameters to 
all measurements listed in Table 1 except ficDM^ 2 , for 
which only fi t heoh 2 < ^ob&h 2 is required. The minimum 
X 2 value is 20.4 for 21 degrees of freedom, corresponding 
to aP-value of 49.3%. 



Parameter 


Best Fit 




Uncertainty 


sign(/j) 


+ 1 






a* 


0.1176 


± 


0.0025 


l/oi em 


127.925 


± 


0.020 


m z (GeV) 


91.1866 


± 


0.0021 


mt (GeV) 


172.2 


± 


1.1 


tan p 


9.0 




+ 11.4 

-3.6 


M 1/2 (GeV) 


303.4 




+ 133 
-59.0 


M (GeV) 


27.6 




+ 122 

-27.5 


Ao (GeV) 


143.2 




+850 
-478 



Predicted Mass Spectrum of SUSY Particles LE n, he0 <!2 . 



p 1800 
O 1600 
S 1400 



1200 



| 1a Environment 
2a Environment 
■ Best Fit Value 





h°A H°H + xJz°Z°X°x;x;T R T L : f 1 T 2 q R q L b 1 b 2 f 1 t 2 9 

Fig. 9: SUSY mass spectrum as predicted by mSUGRA 
parameter fit to low energy measurements, requiring 
^theoh 2 < Q bsh 2 with sign(/i) fixed to +1. 



n ohs h 2 and contributes (f2 thco h 2 - S7 ohs h 2 )/a ncDMh 2 for 
fitheoh 2 > Slobsh 2 . The result of this fit is shown in Ta- 
ble 7. The overall x 2 of the fit improves marginally with 
respect to the baseline fit. The preferred regions of tan/3, 
M and A move to slightly smaller values, but the varia- 
tions are not large compared to the uncertainties. AIq re- 
mains almost unchanged. The uncertainties increase with 
respect to the baseline fit, but the order of magnitude re- 
mains the same. The two-dimensional 95 % CL areas of the 
parameter projections are increasing, but no qualitatively 
new features are observed, as seen in the overlay of param- 
eter regions with the baseline fit in Figure 4. Therefore, for 
this modification of the baseline fit, the predicted collider 
observables do not differ significantly from the baseline fit. 

In Figure 9, the predicted Higgs boson and sparticlc 
mass ranges are shown. The observed mass ranges are very 
similar to those observed in Figure 6. While the overall 
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Table 8: Result of the fit of the mSUGRA model with 
sign(^) = +1 including four additional SM parameters to 
all measurements listed in Table 1 except J7cdm^ 2 - The 
minimum x 2 value is 20.4 for 21 degrees of freedom, cor- 
responding to a "P-value of 49.3 %. 



Parameter 


Best Fit 




Uncertainty 


signO) 


+ 1 






a s 


0.1177 


± 


0.0019 




127.924 


± 


0.014 


m z (GeV) 


91.1870 


± 


0.0022 


mt (GeV) 


172.2 


± 


0.83 


tan p 


10.9 




+9.9 
-5.1 


Mi /a (GeV) 


316.2 




+ 122.9 
-69.7 


Mo (GeV) 


45.1 




+ 119.3 
-43.8 


A (GeV) 


209.1 




+973.5 
-494.1 



Predicted Mass Spectrum of SUSY Particles LE no £1 | 
2000- 



> 1800 
£ 1600 
8 1400 



■i= 1000 



| 1a Environment 
2a Environment 
' Best Fit Value 




h° a°h° h + x ° x° x° %\ x; r 2 T R T L h h % q L 6, b 2 t, f 2 g 



Fig. 10: SUSY mass spectrum as predicted by mSUGRA 
parameter fit to low energy measurements except ^?cdm^ 2 
with sign(/j,) fixed to +1. 



range of masses remains similar, there are subtle changes 
in the mass differences, which are explained below. 

Removing i7cDM^ 2 completely from the list of observ- 
ablcs does not change the fit result significantly. The best 
fit point is identical within the 1 a uncertainty with the 
best fit point for i7 t hco^ 2 < ^obsh 2 . The fit result is pre- 
sented in Tabic 8. Again the uncertainties increase with re- 
spect to the baseline fit and the fit with i7 t heo^ 2 < ^obsh 2 , 
but the order of magnitude of the results and single pa- 
rameter uncertainties remain unchanged. This is exempli- 
fied in Figure 11, which shows the direct comparison of 
the fit results with and without ilcDMh 2 . The increase in 
correlation through the addition of J?cdm^ 2 is clearly vis- 
ible. The predicted mass ranges in Figure 10 are similar 
to Figure 9. 

While the overall range of the accessible parameter 
space is not affected very strongly by including or exclud- 
ing the constraint on I2cdm^ 2 in the presence of the other 



2D 95% CL no n, all LSP 
1D68%CL no a, all LSP 
2D95%CL no n, LSP y° 
1D68%CL non, LSPx 4 
2D 95% CL £lth eo <n obs ' 
1 D 68% CL fl, heo <a. bs 
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200 



300 400 
m, (GeV) 



500 
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Fig. 12: Region in the space of the LSP mass m^o and 
rrifj , which for most parameter points is the NSLP, as 
predicted by a mSUGRA model compatible with existing 
data. Using all constraints, only a area with small tonlsp — 
wilsp remains. 



strong obscrvables from Table 1, there is a strong influ- 
ence of J?cdm^ 2 on the mass differences between the spar- 
ticles. This is shown in Figure 12 using the prediction for 
ra^o and 1 . If no requirements on the LSP are made at 
all, i. e. if also charged LSPs are allowed, a large area of 
chargino and stau masses are allowed. The requirement of 
a neutral LSP excludes the upper left area in Figure 12, 
where mf 1 < to^o . The remaining area is allowed with- 
out applying further constraints on /?cdm^ 2 - It exhibits 
a wide range of mass differences between LSP and NLSP, 
between GeV and 450 GeV. The additional constraint 
of /2 t hooft- 2 < ^obsh 2 severely constrains this mass differ- 
ence. Only a small band close to i? t hco^ 2 = i^obsh 2 re- 
mains, with — m^o < 22 GeV at 95 % CL. As already 
outlined in the discussion of the baseline fit, this result 
has significant impact on the expected exclusive collider 
measurements. This shows that the correlations among 
the parameters changes in a significant way if ficuMh 2 is 
included into the fit, while the overall range of parameter 
uncertainties remains similar. 

The result of a mSUGRA fit to the observables from 
Table 1, but with < xp = af™ is shown in Table 9. 
This shows the allowed mSUGRA parameter space for 
sign(/i) = +1 in case the difference between measurement 
and theoretical SM prediction of = (g — 2)^/2 should 
vanish in the future, either due to statistical fluctuations 
or due to systematic shifts in the prediction because of a 
different treatment of the hadronic corrections. The pre- 
ferred parameter space in this case changes dramatically 
with respect to the baseline fit, in contrast to the fits with 
different treatments of i?cDM^ 2 - The best fit value of M1/2 
moves to 903 GeV, but with large uncertainties towards 
smaller values. In contrast, the preferred value of Mq re- 
mains small, but with very large uncertainties towards 
larger values. The parameter tan/3 does not change dra- 
matically, and Aq is basically unconstrained. 
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Table 9: Result of the Fit of the mSUGRA model with 
sign(/i) = +1 including four additional SM parameters to 



,SM 



all measurements listed in Table 1, but with a° xp 
The minimum x 2 value is 19.7 for 22 degrees of freedom, 
corresponding to a "P-value of 60.2%. 



Parameter 


Best Fit 




Uncertainty 


sign(/x) 


+ 1 








0.1174 


± 


0.0021 


1/ctem 


127.927 


± 


0.017 


m z (GeV) 


91.1874 


± 


0.0020 


m t (GeV) 


172.50 


± 


1.19 


tan f3 


7.4 




+9.1 
-4.6 


Mi/a (GeV) 


903.6 




+67.5 
-778.4 


Mo (GeV) 


180.1 




+ 1566.0 
-110.0 


A (GeV) 


956.9 




+1515.6 
-3048.3 



Predicted Mass Spectrum of SUSY Particles SM (g-2) 

4UUU I 




„ 1500 
a 

.2 1000 
■a 
a 

a. 500 



h° A H° H* z ° jc° z ° z ° x ; x 2 T R T L % \ \ \ B 2 t, t 2 § 



Fig. 13: SUSY mass spectrum as predicted by mSUGRA 
parameter fit to low energy measurements with a° xp = 
a^ M and with sign(/i) fixed to +1. 



This has significant consequences for the expected 
sparticle masses, as shown in Figure 13. In comparison 
to the prediction for the baseline fit, the situation changes 
dramatically. The gluino and squarks tend towards masses 
between 500 GeV and 3500 GeV, severely reducing the 
expected production cross-section of colour-charged spar- 
ticles. In contrast to the heavy squarks, the slepton and 
electro-weak gaugino masses remain in the range below 
500 GeV. 

The change with respect to the baseline fit is also vis- 
ible in the allowed parameter space in the detailed two- 
dimensional projections, as shown in Figure 14. Two dis- 
tinct parameter regions both for the 1 a and the 2 a en- 
vironment can be seen in several parameter projections, 
together with a tendency towards larger values of the 
mass parameters M and M 1 / 2 and hence higher values 
of the expected sparticle masses, especially for the case of 
sign(^) = -1. 



If the positive deviation of the measurement of 
from the SM prediction is removed, there is no remaining 
clear preference among the measured observables for ei- 
ther sign(/i) = +1 or — 1. The fit of the mSUGRA scenario 
with a° xp = a^ M and sign(/i) = +1 yields a minimal x 2 of 
19.7, while the fit with sign(,u) = —1 yields Xmin = 19-1 
for 22 degrees of freedom. The results of the latter fit are 
shown in Table 10. The best fit values of M ± / 2 and Mo are 
strongly different between the two results, but the uncer- 
tainties span similar areas. No significant differences are 
obtained for tan/3 and A$. 

The expected sparticle masses and uncertainties for 
the mSUGRA fit with aJfP = a™ and sign(^i) = -1 
are shown in Figure 15. In contrast to the result for 
sign(/i) = +1, very light gauginos are expected due to the 
low value of M . The best fit result for the squark masses 
is at around 1800 GeV at a similar level, but the uncer- 
tainty does not extend to much higher values. Especially 
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Fig. 14: Expected mSUGRA parameter regions for the fit to low energy measurements with a 1 
sign(/i) = +1 and sign(/i) = —1 arc overlaid. The best fit point is shown for the result with sign(/i) = +1 



Results for 



Table 10: Result of the fit of the mSUGRA model with 
sign(/i) = — f including four additional SM parameters to 
all measurements listed in Table f , but with a° xp = a,, M . 
The minimum x 2 value is 19.1 for 22 degrees of freedom, 
corresponding to a V-value of 64.0%. 



Parameter 


Best Fit 




Uncertainty 


sign(/i) 


-1 






Cts 


0.1175 


± 


0.0018 


1/ctem 


127.929 


± 


0.018 


m z (GeV) 


91.1872 


± 


0.0022 


mt (GeV) 


172.25 


± 


1.14 


tan f3 


11.65 




+8.03 
-7.92 


Mi /a (GeV) 


129.9 




+819.1 
-28.8 


M (GeV) 


1760.1 




+24.4 
-1654.7 


A (GeV) 


62.1 




+2000.9 
-2016.3 



Predicted Mass Spectrum of SUSY Particles SM (g-2) signu = 

2400 h 
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Fig. 15: SUSY mass spectrum as predicted by mSUGRA 
parameter fit to low energy measurements with a® xp = 
a^ M and with sign(/i) fixed to —1. 



for the sleptons there are two distinct regions predicted, 
which is a consequence of the prominent distinct allowed 
regions in the Mq and Mi/ 2 distributions shown in Fig- 
ure 14. Again, the gluino is expected to be lighter than 
the squarks, leading to decay signatures strongly different 
from the chain of several two-body decays expected for 
the baseline fit. 

In order to explore the significance of (g — 2) M for the 
fit result, Table 11 shows the result of the mSUGRA fit 
with sign(/x) = +1 to the observables from Table 1 ex- 
cluding (<? — 2) M . Due to the lower minimal \ 2 °f 19-5 



than the fit with 



,SM 



no secondary minima are 



observed, therefore smaller allowed regions in Mo are ob- 



tained for the fit without a° xp than with a° xp = a^ M . The 
comparison of Table 11 with the results of the baseline fit 
in Table 3 shows that (g — 2)^ represents an important 
constraint, but that the remaining observables still favour 
SUSY in the mass range below 1 TeV. 

For illustration of the difference between the SM pre- 
dictions of (g — 2) M from e + e~ and r data, the fit of 
mSUGRA using the predicted SM value and uncertain- 
ties from [94] is shown in Table 12 and the corresponding 
mass spectrum in Figure 16. The mean value of the spar- 
ticle masses and mSUGRA parameters is compatible with 
the result of the baseline fit within their uncertainties, but 
the 95 % CL area of the predicted sparticle masses ranges 
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Table 11: Result of the fit of the mSUGRA model with 
sign(^) = +1 including four additional SM parameters 
to all measurements listed in Table 1 except (g — 2)^. 
The minimum \ 2 value is 19.47 for 21 degrees of freedom, 
corresponding to a "P-value of 55.5%. 



Predicted Mass Spectrum of SUSY Particles LE no (g-2) mSUGRA 



Parameter 


Best Fit 




Uncertainty 


sign(^i) 


+1 






a s 


0.1177 


± 


0.0019 




127.925 


± 


0.015 


m z (GeV) 


91.1875 


± 


0.0020 


mt (GeV) 


172.5 


± 


1.1 


tan f3 


7.5 




+8.9 
-4.6 


Mi /a (GeV) 


389.2 




+568.9 
-117.5 


Mo (GeV) 


72.2 




+ 145.0 
-22.1 


Aq (GeV) 


270.0 




+ 1492.5 
-1985.9 
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Fig. 18: SUSY mass spectrum as predicted by mSUGRA 
parameter fit to low energy measurements except (g— 2)^ 
with sign(/x) fixed to +1. 



Table 12: Result of the fit of the mSUGRA model with 
sign(/i) = +1 including four additional SM parameters to 
all measurements listed in Table 1, using the predicted 
SM- value of (g — 2) M from [94]. The minimum x 2 value 
is 19.6 for 22 degrees of freedom, corresponding to a V- 
valuc of 60.6%. The fit results of the SM parameters are 
consistent with the other fits in this section. 



Parameter Best Fit Uncertainty 

sign(^) +1 

tan (3 7.9 t^ 9 6 

M 1/2 (GeV) 350.3 t^ 06 

Mo (GeV) 65.9 ±^ 3 8 s 

A (GeV) 100.4 t^ 1 




Fig. 16: SUSY mass spectrum as predicted by mSUGRA 
parameter fit to low energy measurements using the pre- 
dicted SM-value of (g — 2) M from [94] with sign(/x) fixed 
to +1. 



up to around 3 TeV, in contrast to the baseline fit wit hits 
upper limit at 1.6 TeV at 95% CL. 

The impact of the two most stringent observablcs 
I?cdm'* 2 and (g — 2)^ is compared in Figure 17. The 
dark blue areas indicate the two-dimensional 95 % CL al- 
lowed parameter region for the fit without ^2cdm^ 2 , show- 
ing that the remaining observables constrain the allowed 
mSUGRA parameter space to regions below 700 GeV both 
in Mq and M 1 / 2 , accesiblc at LHC. The light blue region 
indicates the allowed parameter space for the fit without 
(g — 2)^. It can be seen that the two constraints are com- 
plementary Without (<j — 2)jj significantly larger values of 
Mq and M 1 / 2 are allowed, but there is a narrower allowed 
band at low Mq and Mi/ 2 , caused by the constraint from 
1?cdm'i 2 favouring the co-annihilation region. 

Since the uncertainties of the mSUGRA parameters 
for the fit without (g — 2) M in Figure 17 extend to much 
larger values of the mass parameters, also the maximal 
values of the sparticle masses allowed at 95 % CL are in- 
creasing significantly with respect to the baseline fit or 
the fit without J7cdm^ 2 - This is visible in Figure 18. This 
result shows again, similar to the result for a® xp = a™, 
that (g— 2)^ plays a major role in constraining the allowed 
mSUGRA sparticle masses to values below 1.6 TeV in the 
baseline fit . 

The previous results show that J?cdm^ 2 has a strong 
role in constraining the size of the uncertainties, but does 
not affect the best fit result or the shape of the predicted 
sparticle mass spectrum. On the other hand, (g — 2) M has 
a decisive impact on the shape of the allowed parameter 
regions, the best fit values of the parameters, their uncer- 
tainties and hence on the predicted particle spectrum. The 
prediction of a rich sparticle spectrum in the kinematic 
range accessible by the LHC however remains stable for 
all explored fits including or excluding cither J?cdm^ 2 or 
(g — 2) M , which is a very encouraging result for the dis- 
covery potential of LHC. The question remains whether 
there are other observables among those shown in Table 1 
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Fig. 17: mSUGRA parameter regions compatible with all low energy measurements except (a) (g — 2) M and (b) 
^cdm^- 2 for various parameter combinations and sign(/^) fixed to +1. For both cases two-dimensional 95 % confidence 
level and one-dimensional 68 % confidence regions are shown. From the latter la uncertainties for individual parameters 
can be derived from a projection of the area to the respective axis. 



which constrain the mSUGRA parameter space. This is 
addressed by removing both X?cdm^ 2 and (g — 2) M from 
the fit and looking for remaining constraints. The result 
of this fit with sign(/i) = +1 is shown in Table f3. The 
results show that the uncertainties increase strongly with 
respect to the baseline fit as well as with respect to the fits 
without either ^cdm/i 2 or (g — 2) M . This shows both the 
complementarity of the ficDMh 2 and (g — 2)^ constraints, 
respectively, and the lack of other strongly constraining 
observables in Table 1. Since no significant constraints on 
low-energy SUSY can be placed anymore without these 
observables, the corresponding fit with sign(/i) = — 1 is 
omitted. 

A comparison of the size of the allowed parameter 
space in the Mq and M 1 ^ 2 projections, for the fits with- 
out ficDMh 2 , (g — 2) M and without both is shown in Fig- 
ure 19. For the fit without both observables the Markov 
Chains does not completely explore the two-dimensional 
2 a uncertainty space, hence the uncertainty region is to 
be interpreted as a lower bound on the uncertainties. This 
could be remedied by using significantly larger statistics 
in the Markov Chains. For the result shown here this is 




1000 2000 3000 4000 5000 
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Fig. 19: mSUGRA parameter regions compatible with 
the collider precision observables from Table 1 and with 
the requirement of a neutral LSP, but without any require- 
ment on ficDMh 2 and (g — 2)^, overlaid with the allowed 
parameter regions using only no i?cDM^ 2 or no (g — 2)^. 
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Table 13: Result of the fit of the mSUGRA model with 
sign(^) = +1 including four additional SM parameters to 
all measurements listed in Table 1 except /?cdm^ 2 and 
(g — 2)^. The minimum \ 2 value is 19.5 for 20 degrees 
of freedom, corresponding to a V-value of 49.1 %. Due to 
the very weak remaining constraint on the SUSY param- 
eter space, in this fit a complete coverage of the possible 
two-dimensional 2 a parameter space can not be ensured 
within the available statistics. Hence the uncertainties arc 
to be treated as lower boundaries on the uncertainties. 



Parameter 


Best Fit 




Uncertainty 


sign(^) 


+ 1 








0.1179 


± 


0.0021 


1/ctem 


127.925 


± 


0.015 


m z (GeV) 


91.1876 


± 


0.0022 


mt (GeV) 
tan /3 

Mi /a (GeV) 
Mo (GeV) 
Aq (GeV) 


172.4 
7.4 
432.6 
64.3 
387.8 


± 


1.2 

+47.3 

-4.7 

+4405.6 

-307.1 

+6816.7 

-62.8 

+6341.4 

-10364.9 
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Fig. 20: SUSY mass spectrum as predicted by mSUGRA 
parameter fit to low energy measurements except {g — 2) M 
and ^cdm^ 2 with sign(//) fixed to +1. 



not done, since the results show clearly that SUSY within 
the reach of the LHC or the ILC is not ensured without 
using ftcDMh 2 and (g — 2)^. 

A part of the allowed sparticle mass spectrum is shown 
in Figure 20. The 1 a allowed ranges of all particles apart 
from the lightest Higgs boson, neutralino and chargino 
extend well beyond a mass of 4TeV. Hence no discovery 
at LHC or ILC can be predicted apart from the SM-like 
lightest SUSY Higgs boson, which is difficult to be dis- 
tinguished from the SM Higgs boson at the LHC [95]. 
The allowed mass range of the lightest Higgs boson is 
rrih = 1 14. 0^ J 4 GeV, showing that the precision data still 
has a tendency to push the Higgs mass towards the lowest 



values allowed by the direct searches. At the same time, 
the push is not strong enough to exclude that the light 
Higgs boson mass is close to the theoretically allowed up- 
per limit at around m^ hco < 135 GeV in this scenario. 

No detailed effect of other observablcs is tested due 
to the failure of the remaining observables to constrain 
the mSUGRA parameter space to the region accessible 
with the next generation of collider experiments. How- 
ever, this does not mean that the precision observables 
or the flavour physics data will not have a decisive role 
in helping to understand SUSY once sparticles are dis- 
covered. The reasons for the strong impact of the preci- 
sion data are demonstrated in Section 4.3. If a Higgs bo- 
son is discovered, the precision of its mass measurement 
will be much better than the precision of cr mh ~ ±3 GeV 
obtained from the low-energy, cosmological and precision 
data. Therefore a further improvement on the uncertain- 
ties of the precision observables would be beneficial for an 
ultimate cross-check of the collider data below and above 
the electro- weak scale. 

In addition, Figure 20 shows a very clear impact of the 
precision and flavour physics observables and especially of 
the Higgs searches at LEP in form of an implicit lower 
bound on the sparticle masses above or around the cur- 
rent direct search limits (e. g. [90,96,97,98]). The lightest 
neutralino mass is expected above 45 GeV, the lightest 
chargino mass above 100 GeV, and the fi mass above 
110 GeV. A small region with a chargino NLSP is allowed 
at the 2 a level, with most points featuring a fi NLSP. 



4.1.3 Fits of the GMSB Model 

The mSUGRA scenario studied so far has the strong ad- 
vantage that it solves a very large amount of experimen- 
tal and theoretical challenges of the SM. In addition, it 
is very well studied at colliders, making it an ideal test- 
ing ground. Therefore it is the main SUSY scenario stud- 
ied in this paper. However, it is not the only way how 
SUSY can be broken at the expected Grand Unified The- 
ory scale of around vIgut ~ 10 16 GeV. As an example of 
a different SUSY breaking mechanism, Gauge Mediated 
Supersymmetry Breaking (GMSB) [11] is explored. It has 
the disadvantage that it includes a very light gravitino in 
the range of ttIq w 0(1 — lOMeV) as LSP, which leads 
to hot dark matter. This is difficult to be reconciled with 
structure formation in the Universe [99,100]. Therefore, 
I?cdm'* 2 is not included in the analysis. 

In GMSB, there are four continuous variables: tan/3 = 
vijv\ is the ratio of the Higgs vacuum expectation val- 
ues, A is universal mass scale of SUSY particles at the 
GUT scale, M mess denotes the mass scale of the messen- 
ger gauge particles between the SUSY breaking sector and 
the visible sector, and C grav is the scale of the gravitino 
coupling. In addition, there are two discrete parameters, 
namely sign(/i) and the number of messenger fields N$. As 
before, separate fits are performed for different values of 
the discrete parameters. A selection of the results is shown 
in Table 14. It can be seen that there is a similar sensitivity 
to sign(^) through the positive value of a° xp — a^ M . There 
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Table 14: Result of the fit of the GMSB model with different values of sign(/z) and N$ including four additional SM 
parameters to all measurements listed in Table 1 except SlcBWih 2 . The minimum x 2 value is 19.30 for 21 degrees of 
freedom, corresponding to a V- value of 56.5 %. The uncertainties correspond to the entry for sign(/i) = +1 and N§ = 1. 



Parameter 


Best Fit 


Best Fit 


Best Fit 


Best Fit 


Best Fit 


Uncertainty for 

N s = l,sign(/i) = +1 


sign(/x) 


+ 1 


-1 


+1 


+ 1 


+1 




N 5 


1 


1 


2 


3 


4 




tan f3 


19.2 


19.4 


17.9 


18.3 


18.5 


+ 15.3 
-6.7 


A (GeV) 


87050 


307629 


53284 


40080 


32643 


+31970 
-17151 


M mesa (GeV) 


431752 


334662 


688567 


1.038 x 10 6 


539328 


+ 1.74X10 6 
-352952 


Cgrav 


411.4 


446.1 


885.5 


460.1 


3368.1 


+ 10042.5 
-411.3 


Xmin 


19.3 


31.0 


19.4 


19.5 


19.5 



is, however, no sensitivity to N§, since all performed fits 
with iVg = 1, 2, 3, 4 achieve the same value of Xmin = 19-5, 
corresponding to a P-value of 56.5%. It can be seen that 
no sensitive limit can be placed on M mess and C grav , while 
tan/? and A can be constrained. Since the prediction for 
tan (3 and A is stable for different areas in M moss and C grav , 
the remaining large uncertainty does not affect the con- 
strained regions for tan/3 and A given the existing mea- 
surements. This shows an interesting complementarity to 
the expected LHC measurements, where constraints can- 
not be set on tan/3 and C grav or M moss alone, but on 
Cgrav x M mess [101] through the measurement of sparticle 
masses and gaugino lifetimes. 

The pull of the individual variables with respect to the 
best fit result of the fit with sign(/i) = +1 and N$ = 1 is 
shown in Figure 21. The obtained pattern is very simi- 
lar to the pattern for the mSUGRA scenario, which again 
confirms that the tension of the SM with the electro-weak 
precision observables cannot be remedied by Supersym- 
metry, apart from moving the Higgs boson mass close to 
or above the experimental limit. The prediction of GMSB 
is m h = (113.5 ± 2) GeV, very similar to mSUGRA. 

The allowed parameter range in tan /3 and A is shown 
in Figure 22 for different values of N5. It is interesting to 
observe that, as already visible in Table 14, the predictions 
and hence Xmin remains unchanged for different values of 
N5. There are different preferred parameter regions in A, 
leaving A x N$ approximately unchanged. Intermediate 
values of tan/3 are preferred, but neither large nor high 
values can be excluded at the 2 er level. 

The predicted sparticle spectrum is also insensitive 
to ./V5 itself, as shown in Figures 23 to 26, where an al- 
most perfect agreement between the spectra for different 
values of N 5 is shown. Generally the expected spectrum 
looks similar to the expected mSUGRA spectrum. This is 
a strong hint that given the existing measurements, Su- 
persymmetry generally provides for the prediction of a 
rich observable particle spectrum at the LHC, indepen- 
dent of the SUSY breaking mechanism. In GMSB, there 
is a tendency towards higher values of the squark and 
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Fig. 21: Pull for the low energy observables used in the 
GMSB parameter fit with N5 = 1 and sign(ju) = +1, using 
the observables from Table 1 and the best fit point from 
Table 14. 



gluino masses with respect to mSUGRA, but the differ- 
ence in the predictions is not decisive enough to base a 
distinction between the scenarios on the mass hierarchies 
for scenarios with similar visible parts of the decay chains 
in GMSB and mSUGRA. 
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Fig. 22: GMSB parameter regions with sign(^) = +1 Fig. 25: SUSY mass spectrum as predicted by GMSB pa- 
and different values of N 5 compatible with the existing ramcter fit to low energy measurements except J?cdm^ 2 
data from Table 1 except J?cdm^ 2 - with N 5 fixed to 3 and sign(^) = +1. 




Fig. 23: SUSY mass spectrum as predicted by GMSB pa- 
rameter fit to low energy measurements except f^cDMh 2 
with N$ fixed to 1 and sign(/i) = +1. 




Fig. 24: SUSY mass spectrum as predicted by GMSB pa- 
rameter fit to low energy measurements except X?cdm^ 2 
with N 5 fixed to 2 and sign(/i) = +1. 



Fig. 26: SUSY mass spectrum as predicted by GMSB pa- 
rameter fit to low energy measurements except f^cDMh 2 
with N$ fixed to 4 and sign(/i) = +1. 

4.1.4 Conclusions for SUSY Scenarios at Colliders 

The results show that Supcrsymmetry, broken at the GUT 
scale, offers several possibilities to fit the existing preci- 
sion data. In the mSUGRA breaking scenario, all param- 
eters can be constrained and the existing data from cosmo- 
logical, low-energy, flavour physics and precision collider 
sources clearly prefer parameter ranges which are acces- 
sible at the next generation of collider experiments. It is 
shown that as expected I?cdm^ 2 and (g — 2) M provide for 
the most sensitive constraints among the available mea- 
surements. Even if individual, very sensitive, variables are 
removed, or if deviations between data and the SM pre- 
diction are assigned to unknown systematic uncertainties, 
a clear constraint in the accessible mass regions remains. 
Only voluntarily removing both i?cDM^ 2 and (g — 2) M si- 
multaneously from the list of observables used in the fit, 
removes the experimental constraint to the parameter re- 
gion accessible at the LHC. 
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The best fit parameter spectrum and the uncertain- 
ties from present data clearly prefers a SUSY scenario 
with a rich phenomenology both at the LHC and the 
ILC. At the 1 a level, all Higgs bosons, gauginos (apart 
from the gluino) and all sleptons are expected below 
m < 600 GeV. The squarks and gluinos are expected be- 
low m < 900 GeV. While this provides for relatively early 
discovery at the LHC, the rich expected spectrum with 
many concurrent production and decay modes will con- 
tribute to a challenging reconstruction of the LHC ob- 
servables sensitive to SUSY masses and branching frac- 
tions. In particular, a very small mass difference between 
the neutralino LSP and the NLSP, which is the t\ of 
tonlsp — "ilsp < 22 GeV at the 95 % CL level leads to 
dominating decays of the gauginos into final states with 
t leptons. However, the exact branching fractions cannot 
be predicted with strong precision. 

4.2 Expected LHC Measurements 

Based on the results of the previous section, we now in- 
vestigate the prospects for the determination of SUSY 
parameters from future LHC measurements. Within the 
mSUGRA model, the preferred parameters from existing 
LE measurements and constraints clearly point towards 
rather light sparticle masses. In order to be consistent with 
the measured dark matter relic density, co-annihilation of 
the LSP and the NLSP has to contribute to the dark mat- 
ter annihilation process. For this process to be efficient, 
the mass difference between NLSP and LSP has to be 
rather small. For the best fit point, within mSUGRA, the 
difference m(fi) — m(x?) is only 8 GeV, and the differ- 
ence m(eu) — m(x?) is 22 GeV. No detailed experimen- 
tal studies for LHC prospects are available for this spe- 
cific parameter point. However, detailed studies exist for 
the SPSla parameter point with parameters tan/3 = 10, 
A = -100 GeV, Mi/2 = 250 GeV, M = 100 GeV, 
sign(/x) = +1 [68]. These parameters taken at face value 
lead to a significantly larger dark matter relic density 
due to mass differences m{f\) — m(xi) = 37 GeV and 
m(en) — m(x?) = 47 GeV. Apart from this difference, the 
collider phenomenology of SPSla is very similar to that 
of the best fit point. The smaller mass differences lead to 
softer spectra for the final state leptons, which is a small 
caveat to be kept in mind in the following analysis and 
should serve to trigger more optimisation of soft lepton 
identification within the LHC experiments. We assume in 
the following as a plausible scenario that SUSY is realised 
with electro- weak scale parameters derived from the high- 
scale SPSla parameters and will be discovered by the LHC 
experiments. As input measurements we use the observ- 
ables specified in Section 2.2 for three different integrated 
luminosities. If taken at face value, the lightest Higgs bo- 
son mass in SPSla calculated with SPheno is 109 GeV, 
slightly below the LEP exclusion. Given the theoretical 
uncertainty as well as the strong dependence of rrih on the 
top quark mass we do not consider this as an inevitable 
constraint. Technically, we set the LEP Higgs mass limit 
slightly below 109 GeV for the luminosity scenarios where 



Table 15: Result of the fit of the mSUGRA model to the 
expected LHC observables for 1 fb _1 . 



Parameter 


Best Fit 




Uncertainty 


sign(» 


+ 1 






tan /3 


9.1 


± 


3.7 


A (GeV) 


-131.8 


± 


742.1 


M (GeV) 


100.2 


± 


4.2 


M V a (GeV) 


249.7 


± 


6.7 



no Higgs boson has been found yet at the LHC and we 
assume that a 109 GeV Higgs boson could be discovered 
at the LHC with similar sensitivity as for 115 GeV. 

4.2.1 mSUGRA Fit with fixed sign^) 

The good agreement in collider phenomenology be- 
tween SPSla and the mSUGRA best fit point offers the 
possibility to use the wealth of Monte Carlo studies per- 
formed for this benchmark point to attempt a projection 
of the SUSY model discrimination power and parameter 
constraints to the LHC era. This is done by performing 
fits to Toy data which have been obtained by smearing 
the observable values according to a Gaussian around the 
nominal SPSla values as explained in Section 3. In the 
top four plots in Figure 27 the distributions for the fit- 
ted parameters from these Toy Fits are shown assuming 
fj, > for integrated luminosities of 1 fb _1 , 10 fb _1 and 
300 fb _1 . The fitted central values and their uncertainties 
are obtained from a Gaussian fit to the parameter distri- 
butions. The corresponding values are listed in Table 15, 
16 and 17. The two bottom plots of Figure 27 display, as 
examples, the x 2 distributions from the Toy Fits for the 
smallest and largest considered luminosity. The fact that 
the x 2 distributions are in very good agreement with the 
expectations for the respective degrees of freedom provides 
confidence that the x 2 minimisation algorithm works re- 
liably. 

Already with 1 fb _1 of LHC data it is possible to con- 
strain the scalar mass parameter Mq and the gaugino mass 
parameter M1/2 to the level of a few percent due to the 
already relatively precise measurements of the endpoints 
of the Iq and the Uq invariant mass spectra. For the de- 
termination of M1/2 a lso the ttiti measurement in events 
with qjj — > (7X1 decays is important, tan (3 and Aq are more 
difficult to determinee. For tan /3, approximately 40 % pre- 
cision and for Ao only an order of magnitude estimate is 
obtained. With 1 fb^ 1 the best constraints on these pa- 
rameters come from measurements involving third gener- 
ation particles, in particular from mf b . 

The increased precision on Mq for an integrated lumi- 
nosity of 10 fb _1 mainly comes from more precise mea- 
surements of the endpoints of the £q spectra. For M1/2 
also new sensitive measurements become available with 
a larger data sample, in particular m g — m^o. The mea- 
surement of m™ b still remains important to constrain tan (3 
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Fig. 27: Results of the Toy Fits to the expected LHC "measurements" from Table 2 for different integrated lumi- 
nosities. The upper plots show the expected distribution of best fit points, exhibiting the strong increase in precision 
for £ mt > 10fb _1 . The lower row of plots shows the % 2 distributions for the Toy Fits with integrated luminosities of 
1 fb _1 and 300 ft) -1 , showing very good agreement with the expected \ 2 distributions for the respective degrees of 
freedom. 
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Table 16: Result of the fit of the mSUGRA model to the 
expected LHC observables for 10 fb _1 . 



Parameter 


Best Fit 




Uncertainty 


sign(/x) 


+ 1 






tan p 


10.08 


± 


0.84 


A Q (GeV) 


-98.0 


± 


52.9 


Mo (GeV) 


100.1 


± 


2.1 


Afi /a (GeV) 


250.1 


± 


1.2 



Table 17: Result of the fit of the mSUGRA model to the 
expected LHC observables for 300 fb _1 . 



Parameter 


Best Fit 




Uncertainty 


sign(M) 


+ 1 






tan P 


9.98 


± 


0.35 


A (GeV) 


-100.2 


± 


11.1 


Mo (GeV) 


100.0 


± 


0.39 


M 1/2 (GeV) 


250.0 


± 


0.30 



and Aq. In addition m^ ax from xl decays provides valu- 
able additional information on tan/3 and Aq at 10 fb _1 . 
For tan (3 also the ratio of branching fractions (2) starts 
to contribute. 

With 300 fb- 1 of LHC data it will finally be possible to 
constrain Mq and Mi/ 2 down to (a few) pcrmillc level. The 

driving factor is an increased precision on mj^ gh , mj>° w and 
- in case of M 1 / 2 - also on m^ ax . Similarly the improve- 
ment on Aq can be traced back to better measurements 
of those observables which already provide the best con- 
straints for 10 fb _1 , namely m^ ax from xl decays and mf b . 
tan/3 at 300 fb _1 is mainly controlled by measurements 
of quantity (2) and the lightest Higgs mass n%h- For tan/3 
(Aq) a relative precision of approximately 4 % (11 %) is 
finally achieved from the given list of observables. 

Concerning the most constraining observables men- 
tioned above it should be noted that they might be very 
sensitive to small changes of the input measurements. 
Therefore they ought to be taken with some care and 
should not be generalised without further cross-checks. 

4.2.2 Determination of sign(/i) 

For the LHC fit results described above, we have not yet 
discussed how to find the correct sign of fi. Using the 
technique described in Section 3.2, we also checked how 
well the sign of /i can be determined from LHC data. 
This is done by performing Toy Fits for each sign of [i 
to the identical set of Toy data. Figure 28 shows the \ 2 
correlations obtained from such fits. If we choose the 
value for sign of \i which yields the best x 2 for a given set 
of LHC measurements wc can estimate the probability 
to make the wrong choice by counting the number 
of Toy Fits below the red bisecting line in Figure 28 
and normalise it to the total number of Toy Fits. The 
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Fig. 29: x 2 correlations obtained from fits to the same 
Toy data set for two different interpretations of m^ ax . 

corresponding numbers can be read off the plots. Already 
with 1 fb _1 there is a good chance to extract sign(/^) 
correctly. The probability for a wrong choice is less than 
5 %. Based on 10 fb _1 or more of LHC data, sign(/x) can 
be determined with negligible error probability. 

4.2.3 First Investigation of Chain Ambiguities 

All the above LHC fit results are based on the idealised as- 
sumption that one knows the contributing decay chain for 
a given mass spectrum. In reality this is, of course, not the 
case. Therefore there are in general various possible decays 
which can contribute to a specific mass spectrum. In the 
context of this study we do not yet address this problem 
systematically A fully realistic analysis, taking the full 
combinatorics for such chain ambiguities into account, can 
be performed when data are present. Here we constrain 
ourselves to the case where we check the impact of mis- 
interpreting a single observable to test the methodology. 
To accomplish this we perform fits to the same Toy data 
set for each possible interpretation of a mass spectrum. 
As an example we allow for two different interpretations 
of one particular di-lepton endpoint as either (correct) 
TO^ ax (m^o,m^o,m^) or (wrong) TO^ ax (m^o, m^o, m^) 
for the case of 10 fb _1 of integrated luminosity. 

Figure 29 shows the \ 2 correlations obtained from fits 
to the same Toy data set for these two different interpre- 
tations of m^" ax . If one always chooses the interpretation 
which yields the smallest minimal x 2 the probability to 
make a wrong decision is approximately 23 %. Figure 30 
shows the mSUGRA parameter distributions from Toy 
Fits assuming the correct endpoint assignment (blue), the 
wrong interpretation (red) and the distribution which is 
obtained if the one with the lowest minimal \ 2 is always 
chosen (black). It is apparent that this chain ambiguity 
has some impact on the reconstructed parameters leading 
to a bias on the mean and to systematically larger values 
for the uncertainties on the parameters, but these effects 
are rather small compared to the uncertainty on the 
parameters. While this observation certainly cannot be 
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Fig. 28: x 2 correlations for mSUGRA fits with \i > and /i < to the same Toy data set (left/middle/right for an 
integrated LHC luminosity of 1 fbr 1 / 10 fb -1 / 300 ib ^ 1 )- 
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Fig. 30: Parameter distributions (for 10 fb _1 ) of the mSUGRA model of two different interpretations of the data, 
fitted to the same Toy data set. The parameter distributions assuming the correct endpoint assignment is shown in 
blue, those using the wrong interpretation is indicated in red and the distribution which is obtained if always the 
interpretation with the lowest minimal x 2 is chosen is displayed in black. 
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generalised to arbitrary ambiguities in the decay chains, 
the principal method can always be applied. Depending 
on the result of the x 2 comparison, the ambiguity can 
be cither translated into increased parameter errors 
or certain hypotheses can be discarded if they yield 
significantly worse \ 2 - 

4.3 Role of Low Energy Observables in the LHC Era 

Now we address the question to which extent low en- 
ergy measurements still contribute to the determination 
of SUSY parameters once LHC results become available. 
We perform this study for two different SUSY models, 
namely mSUGRA and MSSM18. The input observables 
for these analyses comprise all the low energy observables 
listed in Table 1 in addition to the LHC observables of 
Table 2. To ensure a consistent set of "measurements" for 
these analyses, nominal SPSla values are used for the low 
energy observables instead of the actually measured val- 
ues. 



4.3.1 mSUGRA Fit 

The fit results of the mSUGRA Toy Fits using low en- 
ergy and LHC observables for the three different lumi- 
nosities are shown in Figure 32. The corresponding fitted 
mSUGRA parameters and the corresponding correlation 
coefficients are summarised in Tables 18, 19, 20, 21, 22 
and 23. As described in Section 4.2 tan/? and Aq are 
rather weakly constrained by LHC measurements alone 
at low luminosity. For these parameters the addition of 
low energy measurements clearly improves the situation. 
At 1 fb -1 , (g — 2) M is an important additional measure- 
ment constraining tan j3 and Aq. For Mo, the cold dark 
matter relic density ^?cdm^ 2 becomes the most sensitive 
observable followed by ni^ gh , the most important LHC 
quantity. The precision on Mi/ 2 is still dominated by the 
LHC "measurements" listed in Section 4.2. Nevertheless 
some improvements are also achieved for this parameter, 
mainly due to (<? — 2)^. 

At 10 fb _1 and above the role of low energy measure- 
ments is largely repressed by LHC observables such that 
the precision on the mSUGRA parameters for increasing 
luminosity asymptotically approaches the precision ob- 
tained from LHC observables alone. 

These results can be used to derive the complete SUSY 
particle mass spectrum assuming the mSUGRA model 
and the best fitting parameters. Figure 33 shows the mass 
spectrum for an integrated LHC luminosity of 1 fb^ 1 as 
obtained from low energy and LHC observables. The re- 
spective mean and most probable values are indicated in 
black and red. la, 2 a and 3 a uncertainties are indicated 
by the blue bands. It should be noted that the masses 
derived in this way are model dependent statements and 
not direct mass measurements. While the masses of the 
light Higgs boson, the light gauginos and the sleptons can 
already be constrained quite well, the masses of the heavy 
Higgs bosons, the heavy gauginos and the squarks are still 



Table 18: Result of the fit of the mSUGRA model to low 
energy and LHC observables for 1 fb -1 . 



Parameter 


Best Fit 




Uncertainty 


sign(p) 


+ 1 






tan /3 


10.2 


± 


2.3 


Ao (GeV) 


-76.3 


± 


184 


M (GeV) 


100.6 


± 


3.4 


M V a (GeV) 


250.2 


± 


5.3 



Table 19: Correlation coefficients for the fitted parame- 
ters of the mSUGRA model to the expected low energy 
and LHC observables for 1 fb -1 . 





tan f3 


A 


Mo 


Mi/2 


tan j3 


1.000 


0.534 


-0.405 


0.793 


A 


0.534 


1.000 


-0.184 


0.493 


M 


-0.405 


-0.184 


1.000 


-0.077 


Mi/2 


0.793 


0.493 


-0.077 


1.000 



quite imprecise. The situation improves significantly if one 
goes to a luminosity of 10 fb -1 (see Figure 34). Increas- 
ing the luminosity to 300 fb -1 (Figure 35) again means a 
clear increase in precision with respect to the 10 fb -1 re- 
sult. This means that stringent spectroscopic tests of the 
mSUGRA model will be possible using the experimentally 
accessible sparticles and precise mass predictions are fea- 
sible for those SUSY particles which cannot be directly 
probed at the LHC. 

In addition to Toy Fits we also perform a Markov 
Chain analysis. Figure 36 shows the quantity Ax 2 = 
— 21n(£) + 21n(£ max ) for all possible mSUGRA param- 
eter pairs for the three considered LHC luminosities of 
1 fb-710 fb-7300 fb" 1 (left/middle/right). C is the two- 
dimensional profile likelihood and £ max the global maxi- 
mum of the likelihood. The black dotted contours repre- 
sent A\ 2 = 1 contours. The results are in good agreement 
with those obtained from the Toy Fits and nicely show the 
partly strong correlation between the parameters which is 
also reflected in Tables 19, 21 and 23. 

For illustrative purposes Figure 37 shows the outcome 
of the same Markov Chain for the parameter pair Aq- 
tan/3 using Baycsian statistics. The lines again indicate 
Ax 2 = — 21n(£) + 21n(£ max ) contours but this time C de- 
notes the marginalised posterior probability (using a flat 
prior probability). Compared to the results derived from 
the profile likelihood the contour lines are more jagged for 
the same Markov Chain length. Apart from these fluctu- 
ations good agreement between the results derived from 
the marginalised posterior probability and those from the 
profile likelihood (shown in Figure 36) is found. 

4.3.2 MSSM18 

So far, we only considered SUSY models with specific 
assumptions on the SUSY breaking mechanism, namely 
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Fig. 31: Discrimination power of the expected low-energy and LHC results from Table 1 and 2 for £ lnt = 1 fb 1 (left) 
and C mt = 10 fb -1 in a 2-dimcnsional plot of the minimal x 2 values of simultaneous Toy Fits of two different models 
(sign(/i) = +1 and sign(/i) = —1). In comparison with Figure 28 a clear increase in separation power for C mt = lfb -1 
is observed. 



Table 20: Result of the fit of the mSUGRA model to low 
energy and LHC observables for 10 fb -1 . 



Derived Mass Spectrum of SUSY Particles mSUGRA LE+LHC 1 fb- 1 



Parameter 



Best Fit 



Uncertainty 



sign(^) 
tan f3 
A (GeV) 
Mo (GeV) 
M1/2 (GeV) 



+1 
10.0 
-99.1 
100.0 
250.1 



± 
± 
± 
± 



0.79 
48.3 
1.9 
1.1 



Table 21: Correlation coefficients for the fitted parame- 
ters of the mSUGRA model to the expected low energy 
and LHC observables for 10 fb -1 . 





tan (3 


^0 


Mo 


M 1/2 


tan (3 


1.000 


0.805 


-0.328 


0.415 


A 


0.805 


1.000 


-0.483 


0.548 


M 


-0.328 


-0.483 


1.000 


0.241 


M 1/2 


0.415 


0.548 


0.241 


1.000 



Table 22: Result of the fit of the mSUGRA model to low 
energy and LHC observables for 300 fb -1 . 



Parameter 



Best Fit 



Uncertainty 



sign(^i) 
tan f3 
A (GeV) 
Mo (GeV) 
M 1/2 (GeV) 



+1 
10.00 
-99.1 
100.00 
250.01 



± 
± 
± 
± 



0.36 
12.0 
0.39 
0.33 
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Fig. 33: SUSY mass spectrum consistent with the exist- 
ing low-energy measurements from Table 1 and the ex- 
pected LHC measurements from Table 2 at C mt = 1 fb -1 
for the mSUGRA model. The uncertainty ranges represent 
model dependent uncertainties of the sparticle masses and 
not direct mass measurements. This is especially visible for 
the heavy Higgs states A, H and H ± , for which no direct 
measurement is expected in the SPSla scenario. 



mSUGRA and GMSB (for the fits to low energy measure- 
ments). As shown in Section 4.3.1 the LHC measurements 
together with LE measurements allow to derive tight con- 
straints on the mSUGRA parameters if sufficient lumi- 
nosity is accumulated at the LHC. In this section we in- 
vestigate if it is possible to relax the strong constraints 
imposed on sparticle masses and couplings by the require- 
ment of a specific breaking scenario. If it will be possible 
to measure the parameters of a more general model, like 
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Fig. 32: Results of the Toy Fits to the low-energy measurements and the expected LHC results from Table 1 and 2 
for different integrated luminosities. The upper plots shows the expected distribution of best fit points, exhibiting the 
strong increase in precision for C mt > 10fb _1 . The lower row of plots shows the x 2 distributions of the Toy Fits with 
integrated luminosities of 1 fb _1 and 300 fb _1 , showing very good agreement with the expected degrees of freedom. In 
comparison with Figure 27, a clear increase in precision is observed, especially for C lnt = 1 fb _1 and for the parameters 
tan/? and M . 
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Fig. 36: A\ 2 = — 21n(£) + 21n(£ max ) contours from Markov Chain for the mSUGRA model using observables from 
Tables 1 and 2. £ is the two-dimensional profile likelihood and £ max the global maximum of the likelihood. The black 
dotted contours represent A\ 2 = 1 contours. The plots are for integrated LHC luminosities of 1 fb -1 /10 fb _1 /300 fb _1 
(left/middle/right). 
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Fig. 37: For comparison with Figure 36, the outcome of the same Markov chain for the parameter pair Ao~ta.n(3 is 
shown using Baycsian statistics. The lines indicate Ay 2 = —2 ln(£) + 2 ln(£ max ) contours with C being the marginalised 
posterior probability (using a flat prior probability). The black dotted contours represent A\ 2 = 1 contours. 
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Fig. 34: SUSY mass spectrum consistent with the exist- 
ing low-energy measurements from Table 1 and the ex- 
pected LHC measurements from Table 2 at C lnt = 10 ft -1 
for the mSUGRA model. The uncertainty ranges represent 
model dependent uncertainties of the sparticle masses and 
not direct mass measurements. With respect to Figure 33, 
a clear increase in precision is observed. 



Fig. 35: SUSY mass spectrum consistent with the ex- 
isting low-energy measurements from Table 1 and the 
expected LHC measurements from Table 2 at £ mt = 
300 ftr 1 for the mSUGRA model. The uncertainty ranges 
represent model dependent uncertainties of the sparticle 
masses and not direct mass measurements. With respect 
to Figure 34, again a clear increase in precision is observed. 



Table 23: Correlation coefficients for the fitted parame- 
ters of the mSUGRA model to the expected low energy 
and LHC observablcs for 300 ft -1 . 





tan P 


A 


Mo 


Mi/2 


tan j3 


1.000 


0.356 


0.178 


0.134 


A 


0.356 


1.000 


-0.266 


0.673 


Mo 


0.178 


-0.266 


1.000 


0.391 


Ml/2 


0.134 


0.673 


0.391 


1.000 



c. g. the MSSM18, at the electro-weak scale, properties of 
SUSY breaking models could be investigated in a bottom- 
up approach [102]. 

The results of a Markov Chain analysis of the MSSM18 
model using LE and LHC observables with an integrated 
luminosity of 300 ft" 1 are shown in Table 25. Most pa- 



rameters can be determined to the level of a few per- 
cent, except for third generation sfermion mass parame- 
ters, the trilincar coupling parameters X T = A T — /i tan /3, 
Xf, = Ab — fi tan /3 and X t = A t — /x cot /3 and the Higgs 
parameters tan/3 and vtla- The precision on the Higgs pa- 
rameters tan/3 and toa suffers from the fact that for the 
analysed benchmark point SPSla the heavy Higgs bosons 
are not directly accessible at the LHC. 

Whereas for the mSUGRA fit to LE and LHC observ- 
ablcs the impact of LE observables is almost negligible for 
300 ft" 1 of LHC data, the situation for the MSSM18 is dif- 
ferent. For some parameters the most stringent constraints 
still come from LE measurements. The most prominent 
examples are B(B — > 57) and (g — 2) M . B(B — ► 57) is 
sensitive to the charged Higgs boson mass m#± which in 
turn is tightly connected to the fitted SUSY parameter 
m A- (.9 — 2) M provides the most sensitive constraints on 
Mf L and X T . 
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Fig. 38: SUSY mass spectrum consistent with the ex- 
isting low-energy measurements from Table 1 and the 
expected LHC measurements from Table 2 at C int = 
300 ftT 1 for the MSSM18 model. The uncertainty ranges 
represent model dependent uncertainties of the sparticlc 
masses and not direct mass measurements. 



The obtained MSSM18 fit result can again be trans- 
lated into a corresponding sparticle mass spectrum. This 
spectrum is presented in Figure 38. Again the masses in 
this Figure are model dependent predictions and do not 
represent direct mass measurements. Compared to the 
corresponding result for the more constrained mSUGRA 
model (Figure 35), some of the sparticle masses have sig- 
nificantly larger uncertainties in the MSSM18. This is par- 
ticularly pronounced for the heavy Higgs boson masses, 
which are - as stated above - not directly accessible at 
the LHC for the considered SUSY benchmark point. 

Although not studied explicitly for the MSSM18, one 
may expect that chain ambiguities may have a larger im- 
pact for this model than in the mSUGRA case. Since 
MSSM18 has more independently adjustable parameters, 
different decay chain interpretations can be more easily 
matched with the model due to the increased flexibility. 



4.4 Low-Energy Observables, LHC and Expectations 
for ILC 

The results of the previous sections show that the expected 
data of the LHC allow to obtain rather precise constraints 
on the mSUGRA parameters once sufficient luminosity is 
accumulated. However, for the MSSM18 scenario, the con- 
straints severely diminish due to the increased theoretical 
freedom. The parameter uncertainties typically increase 
by a factor of 10 or more. Therefore an extrapolation of 
the SUSY parameters from the electro-weak scale to the 
GUT scale is afflicted with large uncertainties, if only low- 
energy, flavour physics, electro- weak precision, cosmolog- 
ical and LHC observables are used. 

The expected measurements at the International Lin- 
ear Collider, however, could dramatically increase the 
experimental precision of the measurements of sparticlc 



Table 24: Result of the fit of the mSUGRA model to the 
existing measurements and to the expected results from 
LHC with £ int = 300 fb" 1 and ILC. 



Parameter 


Nominal value 


Fit 




Uncertainty 


tan f) 


10 


9.999 


± 


0.050 


Mi/a (GeV) 


250 


249.999 


± 


0.076 


M (GeV) 


100 


100.003 


± 


0.064 


A (GeV) 


-100 


-100.0 


± 


2.4 



masses and couplings. In addition to just increasing the 
experimental precision, the ILC is also expected to de- 
liver a wealth of measurements of absolute branching frac- 
tions and cross-sections, many cross-section times branch- 
ing fraction measurements, and many model-independent 
measurements of quantum numbers and CP-properties. 
This expected wealth of data, especially in a SUSY sce- 
nario with a rich phenomenology below a mass scale of 
500 GeV, as predicted by the present measurements in 
Section 4.1, will strongly enhance the knowledge from the 
LHC due to the expected complementarity of ILC and 
LHC results [13]. 

In this section, first the expected precision on the pa- 
rameters of the mSUGRA model is studied, followed by a 
detailed comparison of the results of the MSSM18 fit us- 
ing only LE and LHC data with those obtained using LE, 
LHC and ILC data. Finally, the increase in precision is 
used to predict the cosmic cold dark matter relic density 
ficDivi'* 2 from collider data, from fits excluding i^cDM^ 2 
itself from the list of observables. 



4.4.1 mSUGRA 



Using the same available and expected measurements as 
in the fit using £ int = 300 ftT 1 of LHC luminosity in Sec- 
tion 4.3.1, plus the expected ILC measurements discussed 
in Section 2.3, the fit of the mSUGRA model to the data 
of the SPSla scenario is shown in Table 24. The compari- 
son with the results without ILC in Table 22 shows the in- 
crease in precision by a factor of 5 to 10. However, the pure 
increase in precision for the fit of a high scale scenario is 
not the only improvement using ILC. First, possible devi- 
ations of the SUSY breaking implemented in Nature from 
a given GUT-scale SUSY breaking scenario, involving as- 
sumptions on unification, are much more visible using also 
ILC data. Second, the high accuracy and especially the 
larger variety (covering couplings, mixings, masses, widths 
and quantum numbers) and stronger model independence 
of the measurements allow to fit more general models of 
New Physics. This makes it possible to study the SUSY 
breaking mechanism using a bottom-up instead of a top- 
down approach. 
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4.4.2 MSSM18 

As discussed in Section 4.3.2, the fit of the MSSM param- 
eters at the SUSY breaking scale allows a bottom-up test 
of SUSY breaking and is independent of any assumptions 
about physics at the GUT scale. Section 4.3.2 showed that 
for the MSSM18 model, the parameter uncertainties from 
fits to existing data and expected LHC data are larger by 
at least one order of magnitude with respect to the fits of 
the mSUGRA scenario. 

Table 25 shows a comparison of the parameter uncer- 
tainties of the fits of the MSSM18 model using LE data 
in combination with £ mt = 300 fb -1 of data at the LHC 
(LE+LHC300) and the latter plus the expected ILC re- 
sults (LE+LHC300+ILC). The results of Markov Chain 
Monte Carlo scans and Toy Fits are in good agreement, 
therefore just the Markov Chain result is shown. For most 
parameters, the uncertainties decrease by approximately 
one order of magnitude. Interestingly, the increase in pre- 
cision is not only limited to those parameters which are 
linked directly to observables at tree level. For example 
it is expected that the uncertainties of the gaugino mass 
parameters Mi and Mi are significantly decreased at ILC 
due to the increased precision on the x\/2 and xt masses 
and the additional information from precise measurements 
of cross-sections times branching fractions for different po- 
larisations. Also, the precision of the heavy Higgs sector 
parameter is expected to increase dramatically, since 
the heavy Higgs bosons A, H and H ± are not expected 
to be discovered at the LHC in this scenario, but to be 
precisely measured at the ILC [13]. In contrast to those 
measurements, no additional experimental information is 
obtained on the gluino mass or the heavier squark masses 
at ILC. In any case, with the exception of M^ R , all pa- 
rameter uncertainties improve dramatically. The reason 
for this behaviour is the strong decrease of correlations. 
For example, the &x/2 masses are determined by Mq L and 
, but also by the off-diagonal elements m^X^ with 
Xi) = Ab — fx tan [3. Due to the strong increase in the de- 
termination of [i and tan (3 from the measurements in the 
Higgs sector (where also A^ plays a role in loop effects) 
and the gaugino sector, also the precision of the param- 
eter Mi is strongly improved, although no direct mea- 
surement in the sbottom sector is made at the ILC in 
this scenario. This example highlights the importance of 
precision measurements for the detailed unravelling of the 
SUSY spectrum, and it is an example of the complemen- 
tarity of LHC and ILC. 

The resulting derived spectrum of sparticle masses is 
shown in Figure 39. It represents a very strong improve- 
ment over the results without ILC in Figure 38. The Higgs 
sector exhibits the strongest improvement due to the di- 
rect observation of heavy Higgs states. Apart from the 
squark mass rriq L , solely governed by the parameter Mq R , 
the uncertainties of all other derived masses increase dra- 
matically. 

As a final test of the agreement between cosmology and 
collider data, and as a showcase for the predictive power 
of precision collider measurements, additional fits with- 
out ficDM^ 2 arc performed with and without the use of 
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ILC using the Toy Fit technique. The resulting predicted 
values of J?cdm^ 2 are shown in Figure 40 and compared 
with the present and expected experimental precision of 
fi CDM ft 2 from the WMAP [57] and Planck [103] data. The 
prediction of ^?cdm^ 2 from collider data without ILC in 
the MSSM18 model shows a long non-Gaussian tail down 
to J^cdm^ 2 = 0. The Gaussian core of the distribution is 
one order of magnitude wider than the expected precision 
from the Planck satellite. Therefore witout ILC, the relic 
density constraints inferable from particle physics within 
the MSSM18 model do not match the precision of cosmo- 
logical measurements. 
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Table 25: Results of the Markov Chain MC analysis of the MSSM18 model using low energy observables, expected 
LHC results for £ int = 300 fb" 1 and ILC. 



Parameter 


Nominal value 


ILC Fit 




<3"LE+LHC300 


O LE+LHC300+ILC 


Ms (GeV) 


194.31 


194.315 


± 


6.4 


0.068 


Mi (GeV) 
M? (GeV) 


135.76 


135.758 


± 


10.5 


0.071 


193.52 


193.46 


± 


43.0 


0.33 


Mf (GeV) 


133.43 


133.45 


± 


38.2 


0.35 


M- (Gevi 


527.57 


527.61 


± 


3.4 


0.64 


M,- R (GeV) 


509.14 


509.3 


± 


9.0 


9.0 


M Sr (GeV) 
M- tL (GeV) 


504.01 


504.2 


± 


33.3 


2.4 


481.69 


481.6 


± 


15.5 


1.5 


M tfl (GeV) 
tan/3 


409.12 


409.2 


± 


103.8 


1.6 


10 


10.01 


± 


3.3 


0.29 


H (GeV) 


355.05 


355.02 


± 


6.2 


0.88 


X T (GeV) 


-3799.88 


-3795.1 


± 


3053.5 


46.6 


X t (GeV) 


-526.62 


-526.8 


± 


299.2 


4.7 


X b (GeV) 


-4314.33 


-4252.1 


± 


5393.6 


728.7 


Mi (GeV) 


103.15 


103.154 


± 


3.5 


0.046 


M 2 (GeV) 


192.95 


192.95 


± 


5.5 


0.11 


M 3 (GeV) 


568.87 


568.66 


± 


6.9 


1.65 


m A (GeV) 


359.63 


360.07 


± 


+ 1181 
-99.3 


1.83 



In contrast to that, the result including ILC for 
the MSSM18 scenario achieves a relative precision on 

(0%nM h2 )/(McDM h2 ) of °- 2 % > which is an order of magni- 
tude more precise than the expected Planck accuracy. An 
agreement between the collider result and the cosmolog- 
ical measurement would provide strong hints that SUSY 
LSPs make up the vast majority of dark matter and would 
allow to make predictions for direct dark matter search 
experiments. For comparison, the achievable accuracy on 
the relic density is also shown assuming mSUGRA. The 
uncertainty is improved again by a factor of two. 

In summary, for a SUSY scenario in agreement with 
the present cosmological, low-energy and collider data, the 
ILC would tremendously improve the theoretical under- 
standing of a SUSY model by improving the precision 
of bottom-up determinations of SUSY parameters with- 
out assumptions on unification and breaking mechanisms 
at the GUT scale. The precision would ensure that cos- 
mological implications of New Physics could be predicted 
with a precision significantly better than the current and 
expected cosmological measurements. 



5 Conclusions 

We have performed a comprehensive study of current and 
future uncertainties and correlations of the parameters of 
supcrsymmetric models, i. e. the mSUGRA and GMSB 
model as well as the MSSM18. 

For the case of LE data presently available we confirm 
the results of [44] leading to the conclusion that within 
the mSUGRA model, sparticles are predicted to be light 
enough for an early discovery at the LHC. In particular, 
the squark and gluino masses, which determine the major 



production cross-sections at the LHC are below 1 TeV at 
68% CL and below 1.6 TeV at 95% CL. The most sensi- 
tive measurements are the muon anomalous magnetic mo- 
ment (g — 2)^ and the cold dark matter density J?cdm^ 2 - 
For (g — 2)fj, the results rely on the calculation of the 
hadronic vacuum corrections based on e + e~ cross-section 
data. Sparticlc masses are less constrained for scenarios 
where the SM prediction of (g — 2)^ is closer to its mea- 
sured value. This is currently the case for the prediction 
based on r-data for the hadronic vacuum corrections [94] , 
where the heaviest sparticles are expected below 1.4 TeV 
at 68% CL and 3 TeV at 95% CL. With no deviation of 
(g — 2) M , sparticles are still constrained to lie below ap- 
proximately 2 (3.5) TeV at 68 (95) % CL. A good fit of 
the data (excluding J?cdm^ 2 ) can also be achieved within 
GMSB yielding sparticle masses approximately below 1.2 
(2.0) TeV at 68 (95) % CL. This result shows that the 
feature of light sparticles is not exclusively true within 
mSUGRA, although it may not be true within the gen- 
eral MSSM. Furthermore the LE data and the value of 
ficDM^ 2 in particular point towards a small mass differ- 
ence of the LSP and the NLSP which is fi . The mass of 
the lightest Higgs boson is predicted to be just above the 
exclusion of the LEP experiments. 

For the SPSla parameter point, which provides a 
phenomenology rather similar to the region preferred by 
the LE fit we determined the prospects for parameter 
measurements at the LHC for a complete set of experi- 
mentally accessible and well-studied observables. Within 
mSUGRA, a coarse determination of the parameters can 
already be achieved with an integrated luminosity of 1 
fb _1 . The precision can be significantly improved when LE 
data are combined with the early LHC measurements. For 
300 fb _1 , a precision of better than 1 % can be achieved 
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on Mo and M\/2- The parameter tan/3 (Aq) can be de- 
termined to 3.5 (11) % precision. At high luminosity the 
impact of LE data becomes small. 

For the MSSM18, parameter determination is signifi- 
cantly more difficult and requires larger integrated lumi- 
nosity. Nevertheless, for 300 fb" 1 a decent determination 
of all sparticle masses can be achieved to a few percent 
precision with the exception of t-z, Ti and the heavy Higgs 
bosons. Here the inclusion of LE data still has significant 
impact, in particular in constraining third generation spar- 
ticles. With a linear collider like the ILC operating at up 
to 1 TeV, the MSSM18 can be reconstructed with a preci- 
sion increased by approximately one order of magnitude. 

It should be noted that the bulk region of the MSSM 
as exemplified in the SPSla parameter point is certainly 
favourable for the prospects of parameter measurements 
at both the LHC and the ILC. For the LHC, various 
experimental studies of different parameter points exist, 
however no full analysis of more difficult regions exists 
to date. In particular in regions where long decay chains 
with charged leptons are suppressed, the reconstruction of 
SUSY parameters will be substantially more difficult and 
imprecise. However, given the constraints from LE data, 
such scenarios appear less likely. Given the smaller mass 
difference between sleptons and the lightest neutralino ob- 
served in the fit to the LE data when compared with 
SPSla, a more detailed comprehensive experimental study 
of co-annihilation points at LHC would be beneficial. 

In addition to these quantitative results we proposed 
some new methodological approaches to take care of am- 
biguities in the assignment of experimental observables to 
physical final states. It was shown, that in some cases these 
ambiguities may be translated into the uncertainty on the 
parameters when the ambiguities cannot be resolved sta- 
tistically. We have also shown that the Bayesian and Fre- 
qucntist interpretation of Markov Chain Monte Carlo lead 
to very similar results for fits including LHC data when 
flat priors are used in the Bayesian approach. For fits of 
LE data only, however, the two interpretations do not nec- 
essarily agree. It is observed that the Bayesian approach, 
which includes marginalisation of the hidden parameters 
requires a prohibitive amount of computing power. In such 
cases, only the Frequentist interpretation is exploited. 

In the future, the technologies presented in this paper 
will be applied to a larger variety of models and finally to 
real data from the LHC. In addition, the proposed treat- 
ment of the assignment ambiguities will be extended to 
further possible self-consistent interpretations of the data 
and the resulting effect on parameter uncertainties and 
possible exclusions of assignments will be evaluated. It 
also is expected to be important to evaluate the effect 
of theoretical uncertainties stemming e. g. from missing 
higher order effects and differences between different im- 
plementations of RGE running in more detail. In addi- 
tion to the uncertainties itself, the evaluation of correla- 
tions among theoretical uncertainties could be relevant. 
Finally, if SUSY is realised in Nature, sparticles could be 
discovered before the discovery of a SUSY Higgs boson. 
Therefore, the implementation of present and future lim- 



its on Higgs boson production in arbitrary models of New 
Physics, using [104] could be important. 
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